1 


AD-778  511 

SEMI-MARKOV  MODELS  FOR  EVALUATING 
NAVAL  TACTICAL  DECEPTION  CONCEPTS 

James  M.  Moore 

Stanford  Research  Institute 


Prepared  for: 

Office  of  Naval  Research 


April  1973 


DISTRIBUTED  BY: 


National  Technical  Information  Service 
U.  S.  DEPARTMENT  OF  COMMERCE 

5285  Port  Royal  Road,  Springfield  Va.  22151 


_  -  — — 


_ Unclassified _ 

DOCUMENT  CONTROL  DATA  •  R  &  D  7m  77c p  s/j 

Secuttty  c  iMxstttcMtion  ot  till •,  body  ol  mb t tract  0nd  Indexing  mnnotatton  min.1  be  .nl.red  when  the  uvttmU  ttpotl  l»  clMa*tli*d2 _ 

1  ORIGINATING  ACTIVITY  fCofpOTAI*  BUthof)  li#.  HEPOUT  SECURITY  CLASSIFICATION 


Stanford  Research  Institute 
333  Ravenswood  Avenue 
Menlo  Park,  California  94025 


UNCLASSIFIED 


2t>.  CROUP 


13  REPORT  TITLE 


SEMI -MARKOV  MODELS  FOR  EVALUATING  NAVAL  TACTICAL  DECEPTION  CONCEPTS 


4.  descriptive  NOTES  ( Typ •  of  fport  and  indutiv  data*. 

Technical  Note 

T  au  T ho R(S)  fF/r«f  n«m«, 


James  M.  Moore 


|6  REPORT  DATE 


April  1973 


•«.  CONTRACT  OR  GRANT  NO. 


N00014-71-C-0419 


6.  PROJECT  NO.  RF018-02 


7a.  TOTAL  NO.  OF  PACES  7b.  NO.  OF  REFS 

127  10 

M.  ORIGINATOR’S  REPORT  NUMBER!** 

NWRC  TN-42 


c.  NR  274-008-40 


[tb.  OTHER  REPORT  NOISI  (Any  othat  number*  that  may  ba  aaalanad 
fb/»  raport) 


10  DISTRIBUTION  STATEMENT 


Approved  for  public  release;  distribution  unlimited. 


II  SUPPLEMENTARY  NOTES 


[15  ABSTRACT 


12.  SPONSORING  MILITARY  ACTIVITY 


Naval  Analysis  Programs  (Code  431) 
Office  of  Naval  Research 
Arlington,  VA  22217 


A  semi -Markov  model  representation  of  a  submarine  searching  for  a 
high  value  target  in  a  field  of  decoys  is  extended  beyond  previous  work. 

The  model  has  been  used  to  assess  the  potentials  of  tactical  deception 
concepts  in  modern  antisubmarine  warfare.  The  assessment  results,  including 
model  sensitivity  analyses,  are  published  in  a  separate,  classified,  final 
project  report.  This  technic*;?  note  describes  the  details  of  the  mood's 
structure.  The  advantages  of  the  earlier  model  are  developed,  and  methods 
for  numerical  solution  of  the  newer,  larger  model  are  examined.  Model 
comparisons  and  sensitivity  are  explored.  Suggestions  for  future  extensions 
and  modifications  are  provided. 


!>V 

NATI0NA!  TfCWNiCAl 
information  rfrviai- 


,1473  (PAGE  ;) 


S/N  0101.t07.M01 


Unclassified 

l#curi«y 


Unclassified 


cutity  Claialficatlon 


ROLK  I  »T  I  ROLt  |  *T  I  ROLt 


Search 


Detection 


FORMAC 


Markov 


Semi -Markov 


Decoy 

Submarines 


Antisubmarine 


Tactical  Deception 


«  MOV  «•  1473  (BACK) 

(PAGE  2) 


Unclassified 

Security  Claaiirieation 


Nava /  Warfare  Research  Center 


April  1973 


Technical  Note 
NWRC-TN-42 


SEMI-MARKOV  MODELS  FOR  EVALUATING 
NAVAL  TACTICAL  DECEPTION  CONCEPTS 


By:  JAMES  M.  MOORE 


Prepared  for: 

NAVAL  ANALYSIS  PROGRAMS  (CODE  431) 
OFFICE  OF  NAVAL  RESEARCH 
ARLINGTON,  VIRGINIA  22217 


CONTRACT  N00014-/1-C-0419 
Task  NR274-008-40 


SRI  Project  1318-5 


Approved  by: 

LAWRENCE  J.  LOW,  Director 
Naval  Warfare  Research  Center 


Approved  for  public  release;  distribution  unlimited. 


IH 


ABSTRACT 


A  semi-Markov  model  representation  of  a  submarine  searching  for  a 
high  value  target  in  a  field  of  decoys  is  extended  beyond  previous  work. 
The  model  has  been  used  to  assess  the  potentials  of  tactical  deception 
concepts  in  modern  antisubmarine  warfare.  The  assessment  results,  in¬ 
cluding  model  sensitivity  analyses,  are  published  in  a  separate,  classi¬ 
fied,  final  project  report.  This  technical  note  describes  the  details 
of  the  model’s  structure.  The  advantages  of  the  earlier  model  are  devel¬ 
oped,  and  methods  for  numerical  solution  of  the  newer,  larger  model  are 
examined.  Model  comparisons  and  sensitivity  are  explored.  Suggestions 
for  future  extensions  and  modifications  are  provided. 


PREFACE 


The  research  reported  on  in  this  technical  note  was  conducted  as  a 
subtask  within  a  larger  project  directed  toward  the  assessment  of  tacti¬ 
cal  deception  techniques  in  surface  force  defense  against  cruise-missile 
armed  submarines.  The  project  was  sponsored  by  Naval  Analysis  Programs, 

Mr.  R.  J.  Miller,  Director,  in  the  Office  of  Naval  Research.  Mr.  J.  G. 

Smith  was  the  ONR  Project  Scientific  Officer. 

The  research  effort  was  performed  within  the  Undersea  Warfare  Pro¬ 
gram  (Mr.  A.  Bien,  Manager)  of  the  Naval  Warfare  Research  Center  (Mr.  L.  J. 
Low,  Director)  of  Stanford  Research  Institute.  Mr.  J.  M.  Moore  succeeded 
Mr.  M.  W.  Zumwalt  as  project  leader. 

The  author  wishes  to  acknowledge  the  assistance  of  Mr.  A.  Korsak  in 
resolving  some  computational  difficulties,  and  of  Messrs  S.  Uralli  and 
F.  Rickey  in  implementing  the  FORMAC  computer  language. 


SUMMARY 


Introduction 

Evaluation  of  tactical  deception  concepts  in  modern  naval  warfare 
is  a  fundamental  concern  of  naval  planners.  This  technical  note  presents 
the  current  status  of  an  evolving  family  of  analytical  models  designed  to 
assist  in  these  evaluations.  This  work  supports  a  larger  methodology 
assessment  and  concept  evaluation  presented  in  Reference  1.  The  reader 
should  consult  this  classified  reference  for  results  using  the  models. 

At  the  present  stage  the  models  deal  primarily  with  the  search  stage  of 
encounter.  However,  their  structure  naturally  leads  to  an  expansion  of 
the  engagement  aspects. 

Problem  Definition 

The  scenario  for  analysis  involves  a  high  value  ship  (HVS)  operating 
within  a  certain  specified  area.  Typically,  the  HVS  is  an  aircraft 
carrier  (CV)  operating  alone  or  with  a  carrier  task  group  (CTG) .  Opera¬ 
ting  within  the  same  area  are  a  number  of  low  value  devices  called  decoys 
which  have  characteristics  similar  to  the  HVS  or  task  group.  A  submarine 
enters  the  area  at  time  zero  and  begins  looking  for  the  HVS  in  order  to 
engage  it.  In  the  search  process  the  searcher  encounters  the  decoys. 
Because  the  destruction  of  a  target  will  involve  the  expenditure  of  a 
limited  resource,  as  well  as  the  possible  compromise  of  concealment,  the 
submarine  does  not  wish  to  destroy  decoys;  he  must  therefore  spend  time 
to  classify  a  target  as  high  value  or  decoy.  The  submarine  wants  to 
minimize  search  time,  while  the  HVS  wishes  to  maximize  it.  The  targets 
are  not  able  to  observe  the  searcher. 
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The  following  questions  about  this  situation  are  of  concern: 

(1)  What  risk  of  being  found  does  the  HVS  incur  over  a 
given  period  of  time? 

(2)  What  are  the  relative  effects  on  risk  of  decoy  charac¬ 
teristics,  such  as  number,  detectability,  mobility,  and 
fidelity? 


These  notions  involve  dynamic  probabilistic  systems  and  they  are 
discussed  in  terms  of  semi -Markov  random  processes. 

Approach 

The  problem  of  search  in  the  presence  of  decoys  naturally  lends 
itself  to  representation  by  a  semi-Markov  process  because  the  state  of 
the  searcher  and  the  transitions  between  states  can  be  well  defined. 
Placing  the  problem  in  the  format  of  a  more  general  semi -Markov  process 
provides  a  framework  which  can  accept  hypothetical  or  experimental  data 
that  might  not  be  well  adapted  to  the  more  specialized  Markov  chain.  A 
feature  of  the  approach  developed  for  this  study  is  that  the  solutions 
to  the  various  molels  considered  are  in  closed  form,  amenable  to 
straightforward  calculation.  The  approach  also  provides  a  framework  in 
which  efficient  event-step  simulation  models  can  be  developed.  The 
assumptions  governing  the  structure  of  the  models  are  presented  in  the 
following  pages. 
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ASSUMPTIONS : 

Geometry 

and 

Kinematics 


There  is  a  single  submarine  opposing  several  HVSs 
and  decoys;  typically,  there  is  one  HVS,  which  is 
a  CV  operating  alone  or  in  a  CTG 

All  activity  takes  place  in  a  circular  operating 
area  of  known  fixed  size 

The  submarine  is  initially  in  the  operating  area 
with  no  contacts  at  time  zero 

At  any  fixed  instant  in  time,  positions  of  HVSs 
and  decoys  are  distributed  over  the  operating  area 
according  to  a  uniform  probability  distribution 

Over  short  intervals  of  time,  everything  moves  in 
straight  lines  with  uniformly  random  headings 

Stochastic  independence  of  motion  is  assumed 
among  decoys  as  a  group,  among  HVSs  as  a  group, 
and  between  the  two  groups,  except  in  the  case  of 
reduced  overlap  when  HVSs  avoid  decoys  on  a 
probabiH  stic  basis 


ASSUMPTIONS : 

Target  Detection 
and  Selection 


Submarine  detection  of  an  HVS,  CTG,  CTG  element 
or  decoy  is  represented  by  a  definite  range  law 

The  centers  of  all  detection  circles  lie  inside 
the  operating  area  (i.e.,  the  searcher  as  well  as 
the  HVS  knows  what  the  operating  area  is) 

When  presented  with  an  array  of  targets  (both 
decoys  and  HVSs)  from  which  one  is  to  be  selected 
for  classification,  the  targets  are  all  equally 
likely  to  be  chosen  (i.e.,  the  decoys  are  identi¬ 
cal  and  are  indistinguishable  from  the  identical 
HVSs  until  a  classification  is  made) 


There  are  no  acoustic  masking  effects 


ASSUMPTIONS: 

Target 

Classification 


*  Misclassification— -a  final  declaration  of  the  tar¬ 
get  to  be  other  than  it  really  is — occurs  with  a 
probability  that  may  depend  on  true  target  type 

Once  a  target  is  classified  as  a  decoy,  that  tar¬ 
get  is  removed  from  the  target  field  as  long  as 
information  on  it  is  retained 

Information  on  a  classified  target  is  retained 
until  the  end  of  a  random  length  of  time  or  until 
another  target  is  detected,  whichever  occurs  first 

Once  a  target  is  selected,  the  submarine  stays 
with  that  one  target  until  classification  is  made 

The  classification  process  is  represented  by  a 
random  length  of  time  immediately  preceding  a 
final  decision  on  target  type 

The  final  decision  on  target  type  is  dichotomous — 
HVS  or  decoy 


If  the  engagement  phase  is  considered,  it  is 
assumed  the  submarine  allows  itself  only  one 
attack  in  which  all  its  weapons  are  launched; 

ASSUMPTIONS:  hence,  an  attack  on  a  misclassif ied  decoy  removes 

the  submarine  as  a  threat  to  the  HVS 

Engagement 

*  The  game  ends  when  a  target  (HVS  or  decoy)  is 
classified  as  HVS  (and  the  submarine  launches 
its  weapons)  or  time  runs  out 


I* 


Exponential  holding  times  in  search  states  are  assumed  throughout. 
Holding  times  in  classification  states  may  have  more  general  distribu¬ 
tion.  In  computer  implementation  there  are  actually  two  models — one  has 
misclassif ication  capability,  the  other  has  the  more  general  capture 
time  distribution  capability.  However,  it  is  often  convenient  to  con¬ 
sider  the  two  as  one,  with  the  above  two  capabilities  mutually  exclusive. 

Simulation  studies,  conducted  as  an  adjunct  to  both  the  original 
three-state  model  formulation  and  the  current  study,  have  demonstrated 
that  the  assumed  exponential  search  time  distribution  and  the  adopted 
model  for  determining  mean  time  to  target  detection  are  in  fact  valid 
with  respect  to  the  assumption;)  made  concerning  the  search  process 
(primarily  the  definite  range  detection  law  and  random  motion) . 


Two  computational  approaches  are  considered.  One  is  based  on 
partial  fraction  expansion  of  exponential  transforms  of  probability 
functions.  This  approach  involved  an  experimental  effort  to  use  the 
IBM  computer  language  FORMAC  (p.59)  to  execute  the  algebraic  manipu¬ 
lations  involved.  Partial  fraction  expansion  allows  consideration  of 
non-exponential  holding  times.  When  all  holding  times  are  assumed 
exponential,  a  more  efficient  approach  based  on  a  Pade'  approximation 
for  a  matrix  exponential  is  utilized.  There  are  two  major  reasons  to 
consider  the  partial  fraction  expansion:  (1)  to  study  sensitivity  to 
distributional  form,  especially  variance,  and  (2)  to  study  character¬ 
istics.  other  than  probability  of  state  as  a  function,  such  as  mean  time 
in  state  and  transition  counting. 
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and  capture  tbe  submarine  for  an  average  of  10  hours  per  decoy/ submarine 
encounter.  Probability  of  a  submarine/:iVS  engagement  over  a  seven  day 
period  is  displayed  as  a  function  of  the  number  of  decoys  for  two 
different  decoy  detection  ranges.  With  no  decoys  deployed  the  HVS  has 
probability  one  of  an  encounter  with  the  submarine.  If  five  decoys  are 
deployed,  this  probability  falls  to  0.88  if  it  is  a  small  device  and 
0.51  if  it  is  a  large  device. 


Conclusions  and  Recommendations 

The  models  uiscussed  in  this  document  are  implemented  as  operational 
computer  programs  and  have  been  used  in  concept  evaluation  studies  reported 
in  Reference  1.  A  comprehensive  sensitivity  analysis  of  the  models  is 
presented  in  that  reference.  Experience  in  the  development  and  use  of 
these  models  allows  some  general  conclusions  to  be  drawn  about  the  method¬ 
ology  . 

Two  conclusions  about  the  significance  of  problem  components  are  easy 
to  draw;  (1)  potential  for  target  misclassif ication  is  a  highly  signifi¬ 
cant  part  of  the  problem,  and  (2)  retention  of  information  on  the  last 
target  classified  is  insignificant.  Also  very  significant  is  what  has 
been  discovered  about  the  importance  of  geometry.  Geometry  representation 
can  have  an  appreciable  effect,  but  the  effect  is  due  primarily  to  the 
representation  of  the  dissipation  of  swept  area.  The  following  simplified 
approach  is  suggested  for  future  work: 

(1)  Ignore  boundary  conditions 

(2)  Adopt  simple  modification  of  unbounded  area  detection 
rates  via  approximations  based  on  Figure  3.10  or  an 
extension  of  this  figure. 

A  conclusion  about  method  of  analysis  is  that  the  FORMAC  computer 
language  is  an  interesting  and  powerful  tool  for  analyzing  continuous 
time  semi-Markov  process,  but  it  can  become  cumbersome  and  time  consuming 
in  itself.  Its  use  is  best  restricted  to  a  smaJ-  number  of  states 
(certainly  no  more  than  six),  or  to  cases  of  very  sparse  transition 
matrices,  where  it  can  be  very  useful  in  deriving  compact  closed  form 
solutions.  iJowevs-.  for  a  larger  number  of  states  it  appears  a  discrete 
time  approximation  is  the  better  approach  where  consideration  of  non- 
Markovian  processes  is  desirable.  In  the  work  described  here  it  seems 
a  2-hour  discrete  time  step  would  be  the  minimum  necessary  and  perhaps 
longer,  say  5  hours,  would  be  acceptable. 
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1.  INTRODUCTION 


Evaluation  of  tactical  deception  concepts  in  modern  naval  warfare 
is  a  fundamental  concern  of  naval  planners.  This  technical  note  presents 
the  current  status  of  an  evolving  family  of  analytical  models  designed 
to  assist  in  these  evaluations.  At  the  present  stage  the  models  deal 
primarily  with  the  search  stage  of  encounter.  However,  their  structure 
naturally  leads  to  an  expansion  of  the  engagement  aspects. 

The  scenario  for  analysis  involves  a  high  value  ship  (HVS)  operating 
within  a  certain  specified  area.  Typically,  the  HVS  is  an  aircraft  car¬ 
rier  (CV)  operating  alone  or  with  a  carrier  task  group  (CTG) .  Operating 
within  the  same  area  are  a  number  of  low  value  devices  called  decoys 
which  have  characteristics  similar  to  the  HVS  or  task  group.  A  submarine 
enters  the  area  at  time  zero  and  begins  looking  for  the  HVS  in  order  to 
engage  it.  In  the  search  process  the  searcher  encounters  the  decoys. 
Because  the  destruction  of  a  target  will  involve  the  expenditure  of  a 
limited  resource,  as  well  as  the  possible  compromise  of  concealment, 
the  submarine  does  not  wish  to  destroy  decoys;  he  must  therefore  spend 
time  to  classify  a  target  as  high  value  or  decoy.  The  submarine  wants 
to  minimize  search  time,  while  the  HVS  wishes  to  maximize  it.  The  tar¬ 
gets  are  not  able  to  observe  the  searcher.  The  following  questions  about 
this  situation  are  of  concern: 

(1)  What  risk  of  being  found  does  the  HVS  incur  over  a 
given  period  of  time? 

(2)  What  are  the  relative  effects  on  risk  of  decoy  charac¬ 
teristics,  such  as  number,  detectability,  mobility,  and 
fidelity? 
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These  notions  involve  dynamic  probabilistic  systems.  They  will  be  dis¬ 
cussed  in  terms  of  semi-Markov  random  processes. 

The  models  presented  here  have  been  developed  as  a  subtask  on  a 
larger  project  that  is  directed  toward  evaluating  deception  concepts  in 
surface  force  Jefense  against  cruise-missile  armed  submarines.  The 

* 

principal  document  of  the  larger  project  is  the  final  report  (Ref.  1), 
entitled: 

"Naval  Tactical  Deception  Concepts:  Methodology  Assessment 
and  Concept  Evaluation  (U)." 

Reference  2,  entitled: 

"Simulation  Model  for  Evaluation  of  Naval  Tactical  Deception 
Concepts  (U)," 

and  this  technical  note  are  supporting  documents  for  the  final  report. 
The  work  reported  in  this  document  is  a  continuation  of  earlier  work 
reported  in  Ref.  3,  entitled: 

"Semi-Markov  Models  of  Search  in  the  Presence  of  Decoys." 

Reference  3  concluded  that  the  semi-Markov  process  is  an  effective  frame¬ 
work  in  which  to  study  questions  regarding  the  capabilities  of  decoys 
to  delay  a  searcher  looking  for  a  particular  type  of  target. 

The  problem  of  search  in  the  presence  of  decoys  naturally  lends 
itself  to  representation  by  a  semi-Markov  process  because  the  state  of 
the  searcher  and  the  transitions  between  states  can  be  well  defined. 
Placing  the  problem  in  the  format  of  a  more  general  semi-Markov  process 
provides  a  framewf-k  which  can  accept  hypothetical  or  experimental  data 
that  might  not  be  well  adapted  to  the  more  specialized  Markov  chain. 


References  are  listed  at  the  end  of  this  report. 
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A  feature  of  the  semi-Markov  approach  developed  for  this  study  is  that 


the  solutions  to  the  various  models  considered  are  in  closed  form, 
amenable  to  straightforward  calculation.  The  approach  also  provides  a 
framework  in  which  efficient  event-step  simulation  models  can  be  devel¬ 
oped.  A  discussion  of  general  Markov  and  semi-Markov  processes  is  avail¬ 
able  in  Ref.  4. 

The  previous  work  evolved  a  three-state  model  and  developed  an 
analysis  of  the  value  of  acoustic  countermeasures  in  surface  force  de¬ 
fense  (Ref.  5).  This  three-state  model  is  able  to  provide  some  insight 
into  the  problem,  but  it  has  some  shortcomings,  specifically: 

•  Detection  radii  are  assumed  to  be  small  in  comparison 
to  the  size  of  the  operating  area 

•  There  is  a  requirement  that  either  HVS  speed  or  decoy 
capture  time,  or  both,  be  large  enough  to  dissipate 
the  area  swept  by  the  submarine  prior  to  encountering 
a  decoy 

•  A  classified  decoy  immediately  returns  to  the  target 
field  after  classification,  i.e.,  the  submarine  does 
not  retain  any  information  from  a  classification 

•  Perfectly  reliable  classification  is  assumed 

•  The  implemented  model  assumed  exponentially  distributed 
classification  times,  even  though  the  methodology  of 
semi-Markov  processes  provides  the  mechanism  for  more 
generality. 

The  modeling  effort  reported  on  here  has  concentrated  on  deriving  and 
implementing  solutions  to  these  problem  areas.  The  following  assumptions 
govern  the  structure  of  the  new  model: 


£ 

Time  spent  by  the  submarine  in  classifying  a  decoy. 
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•  There  is  a  single  submarine  opposing  several  HVSs  and 
decoys;  typically,  there  is  one  HVS,  which  is  a  CV 
operating  alone  or  in  a  CTG 

•  All  activity  takes  place  in  a  circular  operating  area 
of  known  fixed  size 

•  The  submarine  is  initially  in  the  operating  area  with 
no  contacts  at  time  zero 

«  At  any  fixed  instant  in  time,  positions  of  HVSs  and 
decoys  are  distributed  over  the  operating  area  ac¬ 
cording  to  a  uniform  probability  distribution 

•  Over  short  intervals  of  time,  everything  moves  in  straight 
lines  with  uniformly  random  headings 

•  Stochastic  independence  of  motion  is  assumed  among  decoys 
as  a  group,  among  HVSs  as  a  group,  and  between  the  two 
groups,  except  in  the  case  of  reduced  overlap  when  HVSs 
avoid  decoys  on  a  probabilistic  basis 

•  Submarine  detection  of  an  HVS,  CTG,  CTG  element  or  decoy 
is  represented  by  a  definite  range  law 

•  The  centers  of  all  detection  circles  lie  inside  the 
operating  area  (i.e.,  the  searcher  as  well  as  the  HVS 
knows  what  the  operating  area  is) 

•  When  presented  with  an  array  of  targets  (both  decoys 
and  HVSs)  from  which  one  is  to  be  selected  for  classi¬ 
fication,  the  targets  are  all  equally  likely  to  be 
chosen  (i.e.,  the  decoys  are  identical  and  are  in¬ 
distinguishable  from  the  identical  HVSs  until  a  classi¬ 
fication  is  made) 

•  There  are  no  acoustic  masking  effects 

•  Once  a  target  is  selected,  the  submarine  stays  with  that 
one  target  until  a  classification  is  made 

•  The  classification  process  is  represented  by  a  random 
length  of  time  immediately  preceding  a  final  decision 
on  target  type 

•  The  final  decision  on  target  type  is  dichotomous — HVS 
or  decoy 

•  If  the  engagement  phase  is  considered,  it  is  assumed  the 
submarine  allows  itself  only  one  attack  in  which  all  its 
weapons  are  launched;  hence,  an  attack  on  a  misclassif ied 
decoy  removes  the  submarine  as  a  threat  to  the  HVS 
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•  Misclassification — a  final  declaration  of  the  target  to 
be  other  than  it  really  is — occurs  with  a  probability 
that  may  depend  on  true  target  type 

•  Once  a  target  is  classified  as  a  decoy,  that  target  is 
removed  from  the  target  field  as  long  as  information  on 
it  is  retained 

•  Information  on  a  classified  target  is  retained  until  the 
end  of  a  random  length  of  time  or  until  another  target  is 
detected,  whichever  occurs  first 

•  The  game  ends  when  a  target  (HVS  or  decoy)  is  classified 
as  HVS  (and  the  submarine  launches  its  weapons)  or  time 
runs  out. 

Exponential  holding  times  in  search  states  are  assumed  throughout. 
Holding  times  in  classification  states  may  have  more  general  distribu¬ 
tion.  In  computer  implementation  there  are  actually  two  models — one 
has  misclassification  capability,  the  other  has  the  more  general  capture 
time  distribution  capability.  However,  it  is  often  convenient  to  con¬ 
sider  the  two  as  one,  with  the  above  two  capabilities  mutually  exclusive. 

Simulation  studies,  conducted  as  an  adjunct  to  both  the  original 
three-state  model  formulation  and  the  current  study,  have  demonstrated 
that  the  assumed  exponential  search  time  distribution  and  the  adopted 
model  for  determining  mean  time  to  target  detection  are  in  fact  valid 
with  respect  to  the  assumptions  made  concerning  the  search  process 
(primarily  the  definite  range  detection  law  and  random  motion) . 
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2.  MODEL  STRUCTURE 


The  search  problem  is  considered  in  terms  of  a  dynamic  system  de¬ 
scribed  by  its  state  variable.  The  system  is  the  submarine,  and  the 
state  variable  of  interest  is  the  occupation  of  its  sensor/ information¬ 
processing  apparatus.  Two  mutually  exclusive  occupations  are  basic  to 
the  analysis:  search  and  classification.  The  further  partitioning  of 
these  states,  and  the  transitions  between  states  is  the  subject  of  the 
following  paragraphs.  First,  the  general  semi -Markov  process  is  re¬ 
viewed,  and  then  the  models  are  developed. 

2.1  The  Continuous  Time  Semi-Markov  Process 

The  semi-Markov  process  (SMP)  is  a  stochastic  process  in  which 
time  is  the  independent  variable,  and  the  dependent  variable  can  assume 
only  a  denumerable  number  of  discrete  values.  A  given  value  of  the 
dependent  variable  is  a  state  of  the  system,  and  the  collection  of  all 
possible  values  is  the  state  space.  Change  from  a  present  state  to  a 
future  state  depends  only  on  the  present  state  and  is  governed  by  the 
transition  probabilities  of  a  Markov  chain,  while  the  time  required  for 
such  a  transition  may,  in  general,  depend  on  both  the  present  state  and 
the  future  state.  Consider  a  space  of  functions  S  =  fS  (t);  t  5  0,  luefl] 

1  (J> 

This  space  is  the  sample  space  and  consists  of  functions  of  time,  with 

the  generic  functional  form  paramete .xzed  by  tneO.  Thus,  a  sample  "point 

consists  of  the  set  of  points  fS  (t);  t  s  0} ;  such  a  sample  point  is 

U) 

called  a  "realization"  of  the  process  and  is  denoted  simply  S  (*).  We 

id 

will  consider  only  a  finite  number  of  states,  with  a  state  denoted  by 
one  of  the  numbers  1,  2,  ...,  N.  Thus,  for  a  fixed  value  of  t,  S  (t) 
is  in  the  set  {l,  2,  ...,  n]  and  specifies  the  state  of  the  system  at 

7  Preceding  page  blank 


that  time.  Hence,  it  is  characteristic  of  the  process  that  each  reali¬ 


zation  S  (•)  is  a  staircase  function;  we  will  assume  each  realization 
ID 


continuous  from  the  right.  A  typical  realization  is  plotted  in  Figure  2.1 
for  a  four-state  SMP.  The  dots  in  the  figure  represent  "events,"  that  is, 
points  whore  a  transition  into  a  state  occurs.  Transitions  from  a  state 
into  the  same  state  may  occur.  If  a  transition  from  state  i  to  state 


j  (i  A  j)  occurs  at  time  t,  then  S  {•)  will  be  continuous  from  the  right 

in 


but  discontinuous  from  the  left  at  t  =  t.  Formally,  we  say  the  process 
is  in  state  j  at  time  t  if  the  last  event  to  occur  was  of  type  j,  where 
Je  tl,  N}. 


S  (t) 
(JU 


FIGURE  2.1  A  TYPICAL  REALIZATION  OF  A  SMP 


A  SMP  may  be  described  in  terms  of  the  transition  probability  matrix 


P  =  [p  ]  of  its  embedded  Markov  chain  and  a  matrix  of  conditional  hold- 
ij 


ing  time  probability  density  functions  H(t)  =  [h  (t)l.  Whenever  the 


process  er.vers  state  i,  the  next  state,  j,  to  which  it  will  move  is 


determined  according  to  state  i's  transition  probabilities  p^,,  ...,  p.„ 

il  IN. 


After  j  has  been  determined,  the  process  holds  for  a  time  t  in  state  i. 

A  J 


The  holding  times  t  are  positive  valued  random  va~*ables,  each  governed 

•J 

by  a  probability  density  function  <pdf)  h  (•). 

“J 

It  is  convenient  to  use  graphical  concepts  in  constructing  and  analyz¬ 


ing  models  of  semi-Markov  processes.  Two  of  these  concepts  are  of  par¬ 
ticular  interest  here.  One  is  the  transition  diagram,  which  is  composed 
of  a  set  of  nodes  numbered  1  to  N  representing  the  states  and  a  set  of 
directional  arcs  linking  node  i  with  node  j  for  each  positive  p  .  An 
example  is  shown  in  Figure  2.2  for  a  three-state  process  with  a  trapping 
state  (a  state  from  which  there  are  no  transitions  ether  than  to  itself). 
Transition  diagrams  are  useful  in  developing  and  describing  the  structure 
of  a  model.  The  other  graphical  concept  of  interest  is  the  flow  graph. 
The  flow  graph  is  composed  of  the  same  nodes  and  arcs  as  the  transition 

possibility  graph,  but  in  addition  each  arc  is  labeled.  The  arc  from 

e  e 

node  i  to  node  j  is  labeled  with  p  h  (s),  where  h  (s)  is  the  exponen- 
tial  transform  of  h  (t),  i.e.. 


Vs> 


00 

■  / 


dt  e  Sth  l  t) 


An  example  flow  graph  is  shown  in  Figure  2.3.  Flow  graph  concepts  are 
useful  in  analysis  of  a  structured  model. 


FIGURE  2.2  TRANSITION  DIAGRAM 
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FIGURE  2.3  PROCESS  FLOW  GRAPH 


Given  the  components  P  and  H(t),  or  the  process  flow  graph,  we  can 
define  the  interval  transition  probabilities  cp..(t).  This  quantity 
cp  .(t)  is  the  conditional  probability  that  the  process  will  occupy  state 

A  J 

j  at  time  t,  given  that  it  entered  state  i  at  time  zero.  Howard  (Ref.  4) 
shows  that 


9.  ,(t)  =  6  w 
ij  ij  : 


i(t>  +  Pik  f  dT  h: 


k=l  o 


where 


for  i, j  =  1, . . .,N  ;  t  2  0 


[1  if  1  =  j 


10  if  i  1  j 


w.(t)  =  1  - 

i 


rv  r  v 

m  ij  4 


(T)dT 
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In  the  modeling  discussed  in  this  document,  only  the  class  of  tran¬ 
sient  semi-Markov  processes  are  of  interest.  A  transient  process  con¬ 
tains  one  or  more  trapping  states.  If  state  N  is  a  trapping  state, 

then  the  interval  transition  probability  cp  (t)  is  also  the  probability 

iN 

the  process  becomes  trapped  in  state  N  at  or  before  time  t,  given  that 
it  started  in  state  i  at  time  zero.  In  this  interpretation  as  a  trapping 


probability,  cp 


ij 


(t)  will  be  our  measure  of  effectiveness. 


2.2  The  Six-State  Model 

The  basic  structure  of  the  six-state  model  is  shown  in  the  transi¬ 
tion  diagram  of  Figure  2.4.  The  states  are: 

State  1.  Entry:  represents  the  initial  search  by  the 
searcher  prior  to  first  encounter 

State  2.  Classify  decoy  (transient):  represents  capture 

of  the  submarine  by  a  decoy  subsequent  to  a  search 

State  3.  Classify  HVS:  self-explanatory 

State  4.  Search  (reduced  field):  repx’esents  search  sub¬ 
sequent  to  classification  of  a  decoy  before  in¬ 
formation  retained  by  the  submarine  on  that  decoy 
has  dissipated 

State  5.  Search  (full  field):  represents  search  with  no 
information  on  decoys 

State  6.  Classify  decoy  (steady  state):  represents  capture 
by  a  decoy  immediately  subsequent  to  capture  by 
another  decoy. 

The  entry  state  (State  1)  provides  a  framework  for  analysis  of 
searcher  entry  into  the  operating  area.  However,  this  particular  aspect 
of  the  analysis  has  not  been  carried  beyond  the  initial  exploratory  phase, 
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FIGURE  2.4  TRANSITION  DIAGRAM  FOR  SIX  STATE  MODEL 


where  it  was  found  the  approach  being  developed  would  be  very  expensive 
computationally  and  have  relatively  small  effect.  The  approach  examined 
is  discussed  briefly  in  Appendix  A.  The  result  is  that  the  searcher  is 
assumed  simply  to  be  in  the  area  initially  without  a  target.  The  entry 
state  is  thus  equivalent  to  the  search-full-field  state  (State  5).  The 
derivation  of  ti:e  holding  time  densities  and  transition  probabilities 
for  these  two  states  proceeds  as  follows.  Define  the  random  variables 


The  size  of  the  state  space  could  easily  be  reduced  from  6  to  5  by 
eliminating  the  redundancy  of  having  two  equivalent  states.  The  re¬ 
duction  is  not  made  here  because  this  technical  note  is  documentation 
for  computer  programs  which  retain  the  two  equivalent  states  as  a  mat¬ 
ter  of  convenience. 
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where  it  was  found  the  approach  being  developed  would  be  very  expensive 

computationally  and  have  relatively  small  effect.  The  approach  examined 

is  discussed  briefly  in  Appendix  A.  The  result  is  that  the  searcher  is 

assumed  simply  to  be  in  the  area  initially  without  a  target.  The  entry 

* 

state  is  thus  equivalent  to  the  search-full-field  state  (State  5).  The 
derivation  of  the  holding  time  densities  and  transition  probabilities 
for  these  two  states  proceeds  as  follows.  Define  the  random  variables 


★ 

The  size  of  the  state  space  could  easily  be  reduced  from  6  to  5  by 
eliminating  the  redundancy  of  having  two  equivalent  states.  The  re¬ 
duction  is  not  made  here  because  this  technical  note  is  documentation 
for  computer  programs  which  retain  the  two  equivalent  states  as  a  mat¬ 
ter  of  convenience. 
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$' r  i kf-WS, 


T  =  time  to  detect  a  decoy 
D 


T  =  time  to  detect  a  HVS, 
H 


where  the  pdf  of  T  is  f  (•)  and  the  pdf  of  T  is  f  (•).  We  are  now  in 
H  H  D  D 

the  context  of  what  Howard  (Ref.  4)  calls  the  competing  process  model. 

The  pdf’s  f  (•)  and  f  (•)  are  the  pdfs  of  the  competing  random  variables. 
H  D 

Using  Howard's  result  and  assuming 


V‘>  *  V 


,  t  2  0  , 


'V 

fQ(t)  =  PDe  ,  t  £  0  , 


it  is  easy  to  show  that 


B  +  B 

HD  PH 


and  p 


B  +  B 

PD  PH 


and  that 


h12(t>  =  h13(t>  "  <SD  +  V' 


+  PH>,; 


t  s  0 


The  interpretation  of  B„  and  Bn  is  that  they  are  rates  at  which  the  HVS 

H  D 

and  decoys,  respectively,  are  encountered.  The  derivation  of  B  and  B_. 

will  be  discussed  subsequently  in  Section  3.3.  Suffice  it  to  say  here 

that  B  =  n  \  ,  where  X  is  the  rate  when  there  is  one  decoy,  and  n  is 
D  D  D  D  D 

the  number  of  decoys;  and  B„  =  n  X  ,  where  \  is  the  rate  when  there  is 

H  H  I!  H 

one  HVS  and  n  is  the  number  of  HVSs 
H 


A  *  Vh 


and  p 
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Since  State  5  is  considered  equivalent  to  State  1  we  have 


52 


=  P 


12 


and 


53 


=  P 


13 


and  also 


V°  *  h53(->  *  V’>  -  V0 


The  search-reduced-field  state  (State  4)  is  similar  to  the  search- 

full-field  state  (State  5)  except  the  number  of  decoys  is  reduced  by 

one  and  an  information  retention  time  T  with  pdf  f  (•)  is  included. 

R  R 

Assuming 


fR(t)  =  XRe 


-X  t 
R 


t  ^  0 


we  can  again  apply  the  competing  process  concept.  Here  \  is  an  input 

-1  R 

and  interpreted  in  the  form  of  \R  as  the  mean  information  retention 
time.  We  obtain 


42 


43 


45 


(n-l)X 

+  n,  X  +  X_ 

D  D 

H  H  R 

n  X 

H  H 

(n-l)X 

+  n  X  +  \ 

D  D 

H  H  R 

X 

R 

(V1)XD  +  Vh  +  XR 


and 


h4i(t)  =  [(V1)xd  +  Vh  +  \h 
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for  n^  £  1  and  i  =  2,  3,  5.  With  n^  =  0  State  4  cannot  be  reached,  so 
transition  probabilities  and  holding  times  are  net  relevant. 

The  "classify  decoy"  states  (2  and  6)  are  complicated  in  terms  of 
transition  probabilities.  The  two  states  differ  in  the  geometrical  con¬ 
siderations  involved  in  specifying  their  transition  probabilities. 

State  6  is  designated  as  "steady  state"  because  it  assumes  the  dissipa¬ 
tion  of  the  area  swept  by  the  submarine  (hence  clear  of  decoys  and  HVS) 
prior  to  encountering  a  sequence  of  decoys.  Figures  2.5  and  2.6  illus¬ 
trate  the  situation.  Assume  for  convenience  of  discussion  that  decoy 
and  HVS  have  the  same  speed  and  detection  range.  In  Figure  2.5  the  sub¬ 
marine  has  encountered  a  decoy  after  a  period  of  searching  with  no  con¬ 
tacts.  At  the  moment  of  encounter  there  is  a  "swept"  circle  (shaded), 


O  Dataciod  Targat 
It  A  Othar  Targati 

FIGURE  2.5  GEOMETRY  AT  TIME  OF  DETECTION 

centered  at  the  position  of  the  submarine,  which  is  free  of  targets 
except  for  the  one  just  encountered  on  the  boundary.  That  is,  the  prob¬ 
ability  of  a  target  in  the  shaded  circle  is  zero.  Actually,  the  circle 
is  an  approximation  of  the  true  swept  area,  which  trails  behind  the 
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submarine  like  a  giant  teardrop.  The  submarine  does  not  know  the  posi¬ 
tion  of  targets  outside  this  circle.  Suppose  the  classification  takes 

time  t  .  Consider  a  hypothetical  target  just  on  the  boundary  of  the 
o 

shaded  circle  (see  Figure  2.6).  Moving  at  speed  v  toward  the  center 

for  a  time  t  ,  the  hypothetical  target  reduces  the  radius  of  the  shaded 
o 

area  by  a  distance  v  t  . 

H  o 
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FIGURE  2.6  GEOMETRY  AT  TIME  OF  CLASSIFICATION 


Hence,  if  v^t^  is  small  enough,  when  the  submarine  is  through  classi¬ 
fying  the  decoy  there  will  be  an  area  (the  reduced  shaded  circle)  where 
there  is  zero  probability  of  targets  existing.  This  is  important  for 
the  following  reason.  Ws  need  the  instantaneous  probabilities  and 

P  where 
H 


PD  =  Pr  [any  given  decoy  is  present) 
P  =  Pr  [any  given  HVS  is  present), 
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where  "present"  means  within  detection  range  of  the  submarine.  The 


basic  idea  of  the  probabilities  as  obtained  for  State  2  is 


P  «3  - 

D  A  -  A. 


p  *,  - 

H  A  -  A„ 


where 


A  =  operating  area 
o 


A  =  area  enclosed  by  the  intersection  of  A  and  the  decoy 
D  o 

detection  circle 


A  .  =  area  enclosed  by  the  intersection  of  A  and  the  HVS 
H  o 

detection  circle 


A^  =  area  of  reduced  shaded  circle  for  decoy 


area  of  reduced  shaded  circle  for  HVS. 


In  State  4  it  is  assumed  that  the  swept  (shaded)  areas  are  dissipated 


completely,  and  hence  A  =  A  =  0.  The  inclusion  of  this  swept  area 

X  « 


in  State  2  and  its  dissipation  in  State  6  has  found  to  be  of  high  sig¬ 


nificance  in  some  important  cases.  The  actual  definitions  of  P  and  P 

D  H 


are  more  complicated  than  those  given  here  and  are  discussed  in  detail 
in  Section  3.1,  along  with  the  detailed  derivations  of  the  transition 


probabilities  pgi  and  p  .  These  details  are  postponed  so  that  overall 


model  structure  can  be  seen  more  clearly. 


The  holding  time  pdfs  for  States  2  and  4  are  arbitrarily  assumed 
to  be  of  the  form 


h21(t)  =  h61(t> 


n-1 


-*v 


(n-1) 


t  ^  0  (i  =  3,4,6) 
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where  p^  5  0  and  n,  a  positive  integer,  are  parameters  of  the  distri¬ 

bution.  This  form  is  known  as  the  Erlang  distribution  and  its  charac¬ 
teristics  are  discussed  in  Appendix  B.  Another  way  of  stating  this 
assumption  is  to  say  the  capture  time  (call  it  T  )  of  a  decoy  has  the 
specified  form.  Capture  time  of  a  target  is  the  time  spent  by  the  sub¬ 
marine  in  classifying  that  target  as  a  decoy  or  HVS.  The  assumption 
about  T  is  made  for  two  reasons:  (1)  at  the  time  of  model  construe- 
tion  no  satisfactory  model  of  the  classification  process  was  available, 
and  (2)  this  form  is  particularly  amenable  to  analysis  and  obtains  a 
variety  of  shapes  via  proper  parameter  selection.  Figure  B.l  in  Appen¬ 
dix  B  illustrates  the  range  of  shapes  that  can  be  achieved  by  changing 

n  for  a  fixed  value  of  p  .  The  mean  of  the  distribution  is  l  1  irres- 

D  D 

pective  of  the  value  of  n.  Hence,  p^  is  the  mean  capture  time  for  a 
decoy. 

State  3,  "classify  HVS,"  is  the  final  state  to  be  explained.  This 
is  the  trapping  state  in  the  six-state  model.  The  only  nonzero  transi¬ 
tion  probability  is  p  =  1.  Because  the  analytic  measure  we  are  using 

Oo 

is  basically  time  of  first  arrival  in  this  state,  it  turns  out  that  the 
form  of  the  holding  time  pdf  h  (•)  is  irrelevant. 


2.3  The  Ten-State  Model 

The  basic  structure  of  the  ten-state  model  is  shown  in  the  transi¬ 
tion  diagram  of  Figure  2.7.  The  basic  difference  between  this  and  the 
six-state  model  is  the  inclusion  here  of  the  possibility  of  the  submarine's 
misclassifying  a  target.  This  possibility  generates  a  requirement  for 
more  states.  The  states  are: 

State  1.  Entry:  represents  the  initial  search  by  the 
searcher  prior  to  first  encounter 

State  2.  Classify  decoy  (transient):  represents  capture  of 
the  submarine  by  decoy  subsequent  to  search 


State  3 *  Classify  HVS:  represents  classification  of  an  HVS 
subsequent  to  search 

State  4.  Classify  decoy  (steady  state):  represents  capture  of 
the  submarine  by  a  decoy  immediately  subsequent  to 
correct  classification  of  another  decoy  or  misclassi- 
fication  of  an  HVS 

State  5.  Classify  HVS  (steady  state):  represents  classifica¬ 
tion  activity  against  an  HVS  immediately  subsequent 
to  the  correct  classification  of  a  decoy  or  mis- 
classification  of  another  HVS 

State  6.  Search  (decoy  reduced  field):  represents  search  sub¬ 
sequent  to  (correct)  classification  of  a  decoy  before 
information  retained  by  the  submarine  on  that  decoy 
has  dissipated  (the  decoy  field  is  reduced  by  one) 

State  7.  Search  (HVS  reduced  field):  represents  search  sub¬ 
sequent  to  misclassification  of  an  HVS  before  infor¬ 
mation  retained  by  the  submarine  on  that  HVS  has 
dissipated 

State  8.  Search  (full  field):  represents  search  with  no  in¬ 
formation  on  decoys 

State  9.  Terminate  (misclassification):  the  game  ends  with 
a  decoy  incorrectly  classified  as  an  HVS 

State  10.  Terminate  (correct  classification):  the  game  ends 
with  an  HVS  correctly  classified. 

Note  that  the  sequence  of  events  modeled  ends  when  the  submarine  classi¬ 
fies  something  as  an  HVS.  This  situation  leaves  something  to  be  desired 
from  the  point  of  view  of  interpretation  since  subsequent  action  may  be 
relevant.  rnvo  approaches  are  available.  It  can  be  considered  that  the 
framework  is  available  on  which  to  build  models  of  action  subsequent  to 
classification.  Or  it  can  be  assumed  the  submarine  launches  all  of  its 
missiles  upon  classification  of  a  target  as  an  HVS.  In  -he  latter  case 
either  the  HVS  is  attacked  or  the  submarine  ceases  to  be  a  threat. 

Most  of  the  comments  made  on  the  six-state  model  apply  directly 
here  also.  State  1  (entry)  and  its  relations  to  State  2  and  3  are 
identical  to  the  other  model.  Specifically, 
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nX.  n  X 

*  D  D  J  *  H  H 

?  12  =  VD  +  Vh  “  Pl3  =  V  +  Vh 


and 


-(n  \  +  n  V  )t 

ht2(t)  =  ht3(t)  =  <Vd  +  VH)e  ’  t  *  0  , 


where  the  superscript  star  (*)  is  used  to  distinguish  ten-state  entities 

from  their  six-state  counterparts.  The  quantities  \  and  X  are  as 

D  H 

discussed  for  the  six-state  model.  State  2,  "classify  decoy  (transient)," 
and  State  4,  "classify  decoy  (steady  state),"  are  very  similar  to  States 
2  and  6  of  the  six-state  model.  In  particular,  the  holding  time  pdf's 
are 


h*l(t)  =  hj  <o  =  ,v'v 


t  ^  0 


(i  =  4, 5,6, 9) 


which  is  the  form  of  h  (•)  and  h  (•)  with  n  =  1.  The  geometry  in  the 
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ten-state  model  is  identical  to  that  of  the  six-state  model,  so  the 
transition  probabilities  p*^  and  p*^  are  adapted  from  p  and  with 

the  required  modification  to  allow  for  misclassification.  Let  p^  be  the 
probability  the  submarine  misclassif ies  a  decoy  it  investigates.  Then 


p*  =  p  ( 1  -  p  )  and  p* 
p2i  P2iv  PMD  P4i 


*  V1  •  W 


(i  =  4,5,6) 


and 


p*  =  p*  =  p 

*29  *4S  *MD 


State  3,  "classify  HVS  (transient),"  and  State  5,  "classify  HVS 
(steady  state),"  are  conceptually  equivalent  to  States  2  and  4,  respec¬ 
tively.  A  different  mean  classification  time  T  = 


CH 


H 


is  allowed  for 


HVSs  so  that 
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took t 


S'2. 


h3i(t> 


=  h!.(t)  = 


51 


=  11  e 


V 


t  S  0 


(i  =  4,5,7,10) 


Letting  p  be  the  probability  of  misclassifying  a  HVS>  we  obta_n  p 

MH  31 


and  p*  in  a  way  that  is  strictly  analogous  to  obtaining  p*  and  p* 


5i 


2i 


li 


Quantities  py  and  p'  are  computed  in  a  manner  similar  to  that  for  p 

3i  5i  2i 


and  with  exceptions  nored  in  Section  3.2.  Then 
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3i 


=  p^a 


3i 


V 


and 


5i 


=  pL<i  -  P.,.) 


5i 


Mr! 


(i  =  4,5,7) 


and 
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Pt 
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=  P. 
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States  6,  7,  and  8  are  search  states.  State  8,  "search  (full  field)," 
is  equivalent  to  State  1  (entry)  for  the  same  reasons  States  1  and  5  in 
the  six-state  model  are  equivalent,  and  the  comments  there  apply  here 
also.  Specifically, 


.  „  ,  'Vo +  -nV* 

h8i(t)  *  Vd  *  W  •  4  =  0 


D  D 


and 


V 


82  n  +  n  X 
D  D  H  H 


and 


n  \ 

H  H 


83 


Vu  *  Vh 


As  before,  is  the  number  of  decoys,  the  single  decoy  encounter 


rate,^  \  the  single  HVS  encounter  rate, t  and  n  the  number  of  HVSs. 
H  H 


States  6  and  7  are  reduced  field  search  states.  State  6,  "search  (decoy 


reduced  field),"  is  the  seme  as  State  4  in  the  six-state  model.  Hence, 


Discussed  in  Section  3.3, 
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,  t  ^  0 


for  n  >  1  and  i  =  2,  3,  8.  Recall  X.  X  is  the  mean  information  retention 


time.  Similarly, 


R 


3*  = 

62 


63 


p* 

68 


(VUXD 


<n  -lH 

+  n  X  +  X 

D  D 

H  H  R 

n  X 

H  H 

(n  -1)X 

+  n  X.  +  X 

D  D 

H  H  R 

X 

R 

<V1UD  +  V»  -  XR 


State  7,  "search  (NVS  reduced  field),"  is  analogous  to  State  6, 
Instead  of  nD,  njf  is  reduced  by  one,  so  that 


h*i(t>  =  [yD  +  +  *„]• 


for  n  a  1  and  i  =  2,  2,  8  and 
H 


•t“»Vfv1)W' 


t  a  o 


n  X 
D  D 


72  n  X  +  (n  -1)X  +  X 

T>  D  H  H  R 


<V1,XH 
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nv\  +  (nH_1>XH  *  XR 
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Here  the  phrase  just  classified  1st  decoy  refers  to  the  first  decoy 
encountered  after  a  search  period  with  no  contacts  and  not  necessarily 
to  the  first  decoy  encountered  during  the  game.  Several  "1st  decoys" 
may  occur  during  any  given  time  period.  Similarly,  "2nd  or  later  decoy" 
refers  to  any  decoy  in  a  string  initiated  by  a  1st  decoy,  where  a  "string’ 
means  the  submarine  goes  immediately  from  one  contact  to  another  without 
an  interval  of  "no  contact." 

In  computing  these  transition  probabilities  we  use  somewhat  the 
same  approach  as  used  in  Ref.  3  to  compute  the  (y^’s.  However,  one  of 
the  crucial  independence  assumptions  used  there  is  shown  here  to  be  un¬ 
necessary. 

Consider  an  area  of  A  square  nautical  miles  containing  n  decoys, 

o  D 

n  HVSs,  and  one  searcher.  Assume  that  at  any  given  instant  the  decoys 
H 

and  HVSs  are  uniformly  distributed  at  random  over  the  area.  The  com¬ 
pletion  of  classification  of  a  decoy  anticipates  an  event.  The  next 
event  to  occur  can  be  either  start  of  search,  start  of  classification  of 
(another)  decoy,  or  start  of  classification  of  an  HVS.  In  order  for  the 

next  event  to  be  start  of  classification  of  an  HVS,  at  least  one  HVS 
* 

must  be  present  and,  if  one  or  more  decoys  are  also  present,  an  HVS  must 
be  chosen  from  the  array  of  targets.  In  particular 

p  =  P(classify  an  HVs|just  classified  1st  decoy) 

=  P(at  least  one  HVS  present  and  pick  HVS  from  field 
of  available  targets[just  classified  1st  decoy) 


"Present"  means  within  detection  range  of  submarine. 
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P(k  new  decoys  present,  i  HVSs  present,  pick 


V1  "h 

-23=EE  , 

k_o  £_■£  HVS|just  classified  1st  decoy) 


n  -1  n 
D  H 


=  \  ^  P(pick  HVsIk  new  decoys,  1  HVSs,  just  classified 

ui  ut  decoy) 

P(l  HVs|k  new  decoys,  just  finished  1st  decoy) 
P(k  new  decoys] just  finished  1st  decoy) 

As  in  Ref.  3  we  assume 

P(pick  HVS|k  new  decoys,  l  HVSs,  just  classified  1st  decoy) 

=  Pr(pick  HVsIk  new  decoys,  l  HVSs) 

_  1 
SL  +  k 

and  also  that 

P(k  new  decoys] just  classified  1st  decoy) 

=  P(k  new  decoys  within  radius  of  submarine] just 
classified  1st  decoy) 


m 


i  n  -1-k 

Pk  (1  -  P  )  ° 

D2  D2 


where  P^  =  P  (any  one  decoy  is  within  radius  R^  of  submarine! just  classi¬ 
fied  1st  decoy).  The  binomial  form  assumes  stochastic  independence  of 

2 

the  decoys.  In  previous  work  P  =  ttR_/A  was  used.  A  better  approach 

D  Do 

is  developed  as  follows. 

P  „  is  given  in  definitional  form  by 
D2  J 


D2 


^  »(rt)  W  -  Vy^l 

L—J  n  (r, )  A 

j=l  <pv  i  o 
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whore 


i  =  range  index 


j  =  angle  index 


r^  =  a  range  from  center  of  operating  area 
cp  =  an  orientation  angle 

J 

n  =  number  of  range  increments 
r 

n  (r  )  =  number  of  angle  increments  at  range  r  . 
cp  i  i 

This  formulation  is  derived  in  the  paragraphs  below.  Let 

E  =  the  event:  any  one  decoy  is  within  radius  RD 
of  submarine 


=  the  event:  submarine  just  classified  1st  decoy 


the  event:  submarine  is  located  distance  r^  from 
center  of  operating  area 


the  event:  submarine  position  at  start  of  classifica¬ 
tion  is  located  at  angle  cpj  from  submarine  position 
at  end  of  classification  (see  Figure  3.1;  measured 
counterclockwise  from  radius) . 


T 

V  *  L 


i  J 


from  the  law  of  total  probability.  Further, 


P(E,ri,cpjlCi)  =  PtEl^.qyCj)  *  P(*Pj Irj'Cj)  •  PtrJ^) 

To  obtain  Pfr^Jc^)  we  assume  the  submarine  position  is  uniformly  distri¬ 
buted  throughout  A  ,  given  that  C  occurs.  It  is  shown  in  Appendix  C 
o  1 

that  the  distance  r  from  the  center  of  A  has  density  function 
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GEOMETRY  AT  END  OF  DECOY  CLASSIFICATION 


k(r) 


=  —  r  0  £  r  £  R 

A  o 

o 

We  consider  a  discrete  approximation  in  20-nmi  increments  starting  with 
distance  10  nmi,  so  that 

w(r  )  =  IOOtt/A 

1  o 

w(r  )  =  [400  +  800(i-l)]n/A  i  =  2 . N 

i  o  R 

where  Nr  =  integer  part  of 

R^/20  +  0.5  and  [r^,  rg,  rg,  . ..]  =  [10,  30,  50,  ...} 

The  density  and  discrete  approximation  are  illustrated  below  in  Figure  3.2 
for  an  operating  area  of  radius  200  nmi.  Then,  Ptr^Jc^)  =  wCr^). 

The  angle  cp  is  assumed  uniformly  distributed  within  the  operating 
area  and  a  discrete  approximation  taken.  Symmetry  is  taken  into  account 
to  reduce  calculations.  The  number  of  elements  in  the  discrete 


FIGURE  3.2  RANGE  DENSITY  AND  APPROXIMATION 


30 


approximation  depends  on  the  range  r  from  the  center  and  is  denoted 

n  (r  ).  The  angle  increment  is  30°,  with  the  first  angle  at  15°,  giving 
CP  i 

possible  angles  {15°,  45°,  75°,  105°,  135%  165°].  Symmetry  accounts  for 

angles  between  180°  and  360°.  Only  angles  such  that  the  distance 

-1  *  -1 
r  =  uu  is  entirely  inside  A  are  included  in  n  (r. ).  Recall  Lt 
2D  D  o  cpi  D 

is  mean  decoy  capture  (or  classification)  time  and  u  is  submarine  speed. 
Hence,  the  distance  r^  is  a  measure  of  travel  by  the  submarine  while 
captured  by  the  decoy.  So 


=  yv 


It  remains  to  calculate  P(e| r^,  cp  ,  C^).  This  probability  is  taken 


to  be 


P(Ejrifcp  = 


vv  -  vw 


where 


A  =  size  of  operating  area  in  sq  nmi  (this  symbol 
o 

alternately  refers  to  the  geometrical  figure 
as  well  as  the  number  of  sq  nmi) 

BD(r_^)  =  the  intersection  of  A0  and  a  circle  of  radius 
centered  on  a  radial  at  distance  r^  from  the 
center  of  AQ  (the  submarine  is  at  r^) 


Bl^ri*  =  ^ntersecti-on  of  Bo(ri)  and  a  circle  of  radius 
r^p  centered  at  a  distance  r2p  and  angle  cpj  from 


the  center  of  B^(r^), 


'lD  *  [°’  “d  '  Vv  1  ■ 


‘ Inside  Aq  If  r[  +  %  +  cos  £  R*. 
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Here  r  is  the  radius  of  the  (perhaps  only)  partially  dissipated  swept 
ID 

area  present  at  the  initiation  of  1st  decoy  classification.  The  quantity 

li  is  decoy  speed.  The  computation  of  B^r^)  and  B^Cr^,  cp^)  involves 

r  ,  r  .  R  and  R  ,  as  well  as  r,  and  co.5  these  areas,  distances,  and 
ID'’  2D'  o  D'  i  j 

angles  are  shown  in  Figure  3.1. 

It  is  useful  here  to  introduce  the  function  Ap(r,  R^,  Rg),  which 
gives  the  overlap  area  of  two  circles  of  radii  R^  and  Rg  whose  centers 
are  a  distance  r  apart.  See  Figure  3.3.  First  define 


Then 


Ap(r,R1,R2) 


-n  [min(Ri,R2>] 

2 

2  R2 

TTR  -  —  [29  -  sin  20  ] 

2  2  1  1 


+  “  [262  -  sin  202] 


if  r  >  R  +  R 
1  « 

if  r  £  |RX  -  R2| 


if  r  >  jR1  -  R2! 


Thus,  B  (r  )  =  A  (r  ,  R  ,  R  ) .  However,  B  (r  ,  cp  >  is  more  involved, 
Di  FioD  iij 

including  some  numerical  integration.  The  computation  is  given  in  flow 
chart  form  in  Figure  3.4.  This  flow  chart  is  in  generic  form,  i.e.,  it 


describes  a  quantity  useful  for  either  decoy  or  HVS.  When  entered  with 
V  riD  811(1  r2D'  thC  output  is  B1(r1>  ®j)‘ 


FIGURE  3.3  OVERLAP  AREA  Af 


A  comparison  between  this  more  involved  form  of  P ^  with  the  pre- 

2 

vious  approximation  of  nRg/A^  is  shown  in  Figure  3.5.  In  this  figure 

PD2  is  shown  as  a  function  of  Rp  and  assuming  A^  is  200  nmi  in 

radius.  It  is  seen  here  that  the  previous  approximation  for  an  of 

200  nmi  differs  by  a  factor  of  about  20  from  the  new  form  with  mean 

capture  time  of  1  hour  and  by  a  factor  of  2  with  a  mean  capture  time  of 

10  hours.  On  the  other  hand,  for  larger  capture  times  and  <  50  nmi, 

the  previous  approximation  is  good.  Numerical  experiments  have  shown 

that  P  0  for  large  values  of  R  (greater  than  the  radius  of  A  ) . 

Da  D  O 

This  happens  because  of  the  inherent  assumption  of  a  previous  condition 
of  no  contact.  Algebraically  what  happens  as  R^  -»  ®  is  «  r  R^ 
which  implies  B^  -  B^  «  0,  Operationally  this  implies  a  predominance 
of  swept  area  hanging  over  from  the  initial  assumption  of  start  without 
contacts  which  is  unreasonable  for  very  large  R  (relative  to  A  ) . 
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FIGURE  3.5  PROBABILITY  P02  VERSUS  DECOY  DETECTION  RANGE  (RD) 

Having  established  P  we  are  finished  with  the  probability 
P(k  new  decoys  )just  classified  1st  decoy) 

In  order  to  complete  the  derivation  of  p  ,  the  remaining  conditional 

MU 

probability 

Pr(A  HVSs|k  new  decoys,  just  finished  1st  decoy) 


E 


is  now  developed.  We  will  show  this  probability  is  actually  independent 
of  k.  First  consider  k  =  1,  i.e., 


V 


P(1  HVS|l  new  decoy,  just  finished  l9t  decoy) 

It  is  convenient  to  drop  the  "l"  from  "l  HVS”  and  ju3t  conrider  the 

situation  of  a  single  HVS  in  A  ;  this  is  consistent  with  our  later  treat- 

o 

ment  of  several  independent  HVSs.  Assume 
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where  p  is  decoy  classification  rarge.  Each  decoy  has  an  A  area 
D  s> 

associated  with  it.  However,  to  keep  notation  simple  we  avoid  as  far 

as  possible  indexing  A  with  respect  to  decoys.  Note  A„  c  A  each  A  . 

s  H  s  s 

Not  shown  explicitly  in  Figure  3.6  but  part  of  the  analysis  is  area  A^, 

which  is  the  intersection  of  A  (r  )  and  a  circle  of  radius  r  centered 

Hi  1H 

at  a  distance  r  and  angle  <p .  from  the  center  of  A  (r  ),  where 
2H  j  Hi 


r  =  max  [0,  R  -  v  u,  ] 
1H  H  H  n  J 


r  =  uu 
2H  H 


A^  as  the  HVS  analogue  of  the  decoy  and  its  calculation  is  shown  in 
Figure  3.6. 

Now  we  are  interested  in  computing  P(h|  r  .  co  ,  F  ).  For  convenience 

i  j  2 

we  temporarily  drop  r^  and  cp^  and  rewrite  "F2"  as  "2  decoys  eB^, "  so 


P(HlF  )  =  PClfj  2  decoys  eB  ) 

£*  U 


In  order  to  derive  this  quantity  we  have  considered  some  geometrical 
concepts.  We  need  also  the  concepts  related  to  overlap  of  influence 
oetween  HVS  and  decoys.  In  Ref.  3  a  queuing  model  of  these  phenomena 
was  introduced.  Let 

X  =  rate  at  which  HVS  encounters  decoys  (discussed  in 
C  Section  3.3) 

X  =  maximum  acceptable  rate  of  HVS  course  change  for  purpose 
o 

of  decoy  avoidance. 


This  range  plays  a  minor  role  in  the  model;  much  more  important  in  the 
model  context  is  decoy  classification  time. 


For  example,  an  HVS  may  be  subject  to  constraints  such  that  it  would  be 

willing  to  change  course  to  avoid  a  decoy  on  the  average  of  only  once  in 

every  3  hours,  in  which  case  \  =  0.2  hr  \  We  can  now  model  the  situa- 

tion  of  a  single  HVS  as  a  single  server  queue  with  customer  balking, 

that  is,  when  the  customer  arrives  if  the  server  is  busy  the  arriving 

customer  leaves.  In  this  case  the  HVS  is  the  server  and  the  customer 

is  a  decoy.  The  customer  arrival  rate  is  X  ,  and  the  service  rate  is 

c 

Xq.  Customer  arrival  corresponds  to  HVS  contact  of  a  decoy.  Customer 
service  corresponds  to  the  length  of  time  after  a  decoy  avoidance  maneuver 
before  the  HVS  is  willing  to  make  another  such  maneuver.  Such  a  queuing 
system  can  be  represented  by  a  two-state  Markov  process,  where  the  States 
0  or  1  correspond  to  t’  e  number  of  customers  being  served  or,  equivalently, 
whether  or  not  the  server  is  busy.  Let  B  denote  "system  in  State  l" 
and  B  "system  in  State  0."  Assume  the  system  is  in  steady  state.  Then, 
from  Ref.  3  we  have 

X 

c 

X  +  X 

o  c 

X 

o 

X  +  X 

o  c 

P(B)  is  interpreted  as  the  probability  the  HVS  is  unwilling  to  maneuver 
to  avoid  a  decoy  at  any  given  time. 

With  States  B  ar.d  B  defined  we  can  write 

P(h}2  decoys  eB^)  = 

P(h|b  &  2  decoys  eB^)  P(b|2  decoys  eB^) 

+  P(h|b  &  2  decoys  eB^  P(§|2  decoys  eB^) 


P(B) 


P(B) 
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We  have 


P(H }b  &  2  decoys  eBD>  =  0 


because  H  implies  B  since  A,  c  A  .  That  is.  an  HVS  would  be  in  A 

H  s  H 


only  if  it  had  chosen  to  overlap  with  a  decoy  and  hence  was  in  State  "b" 


of  the  overlap  model . 


*  (1)  ( 2 ) 

Examine  Figure  3.7.  The  area  A  is  the  union  of  4  and  A  , 

s  s 


where  the  latter  areas  correspond  to  decoys  and  D^,  respectively. 


Now  write 


P(H|B  &  2  decoys  eB^) 


=  P(H|HVSeA*  &  B  &  2  decoys  eB^) 


Pr(HVSeA*|B  &  2  decoys  eB^) 


+  P<h|hVS$A*  &  B  &  2  decoys  eB  ) 


P(HVS^A*  Jb  &  2  decoys  eB^ 


As  above,  P(H|HVSiA  &  B  2  decoys  eB  )  =  0  since  A  ci  A*#  Further,  assume 

H 


P(HlHVSeA*  &  B  &  2  decoys  eB  )  = 


A  -  A 
H  1 


P(HVSeA*|B  &  2  decoys  eB^  =  — 


A  -  A 
H  1 

P(H|B  &  2  deco  .  :B  )  =  — - - 

Pi 

o 


Here  we  have  assumed  independence  of  HVSeA  and  "h"  from  ”2  decoys  e.1  " 


when  the  former  are  also  conditioned  on  B."  In  addition,  assume 


P(Bl 2  decoys  eB^)  =  P(B),  so  that  we  finally  have 


FIGURE  3.7  AREAS  IN  QUEUE  MODEL  OF  INFLUENCE  OVERLAP 


P(H j  2  decoys  eB  )  =  — - -  P(B) 

D  A 

o 

But  this  indicates  "h"  is  independent  of  "2  decoys  eB^"  regardless  of 
conditioning  on  "B."  Furthermore,  it  is  evident  that  the  same  is  true 
for  "n  decoys  sB^"  all  n  =  l,  2,  so  in  general 


:/^^fr'^V*s  *•*”' VZa</h  '  y  <*^i!|^B 


or  in  terms  of  our  earlier  notation 


P(K|rt»  ,F2)  =  PWj^,^)  = 


AH(ri>  iVvV 


X  +  X 
o  c 


Retreating  back  through  all  the  notation  we  can  now  express 


P(HVs|k  new  decoys,  just  classified  1st  decoy) 

X  «— «  «—«w(r)  A  (r  )  -  A  (r  ,cp  ) 


V  >m  m  m  ■  <»  Wll  /  A  Vf  /  M  U  .V  j 

^ _  y  i  H  i  1  i  Tj 

+  X  Zj  L-J  n  (r. )  A 

c  i  j  cp  i  o 


£  P 


H2 


In  particular,  '  HVS"  is  independent  of  k,  the  number  of  new  decoys. 
Assuming  the  HVSs  are  stochastically  independent,  we  have 


P(j&HVSs|k  new  decoys  &  just  classified  1st  decoy) 

fn.. 


"V  a-?  >v* 

1  H2V  H2 


thus 


.-1 


23 


X  +  k 


V  k  '  Pd2>  (  j-'1 


PH2> 


n  -l 
H 


To  obtain  P ^  we  simply  note 


p„  =  P(resume  search  I  just  classified  1st  decoy) 
24 


=  P(0  HVS  and  0  new  decoys  present] just  classified  1st  decoy) 

=  P{0  HVS  present] just  classified  1st  decoy) 

•  P{0  new  decoys  present | just  classified  1st  decoy) 


- [i  -  -.A1  -  vvl 
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probabilities.  Strictly  speaking,  the  p*  and  p'  are  complete  sets  of 

3i  5i 

transition  probabilities  in  that  E  p'  =  1  and  Z  p'  =1;  they  are  aug- 

i  3i  j[  5i 

mented,  however,  with  p  to  obtain  p*  and  p*  ,  as  described  in  Section  2, 

mh  3i  5i 

The  geometrical  situations  for  p*  and  p*  are  identical  to  those  for 

oi  5i 

P2i  and  P6i>  respectively.  The  difference  in  the  situations  is  that,  for 
the  former,  an  HVS  has  just  been  classified,  while  in  the  latter  a  decoy 
has  just  been  classified.  This  difference  merely  requires  the  substitu¬ 
tions  of  n^  for  n0  ~  1  and  n  -  1  for  n  in  the  computations  for  p 
3nd  P6i  in  order  to  obtain  p^  and  p'  .  Hence,  we  have 


“D  "H*1 


'■■SS*(>-vv 


n  -k/n  -I 


.  n  -l-l 

P*  (1  -  P  )  H 
H2  H2 


n  -1  n 

-  Cl-pJ  t.-PD23D 


p*  =  1  -  p'  -  p' 

*35  *34  *37 


j*D  *  /n\  n  -k  /n  -1\  .  n 

e  Erh  ;k «-vD  : ko-vH 

k=l  1=0  Vk/  \l  f 


n  -1  n 

-  ti  -  V  ti  -  pD6] 


n  -1-& 


=  1  -  p'  -  p' 

*54  *57 


where  P  ,  P  m  P  ,  and  P  „  were  defined  in  the  preceding  subsection. 
D2;  D6  H2*  H6 


3.3  Detection  Rates 

In  developing  the  model  structure,  reference  was  made  to  random 
variables 


T  =  time  to  detect  a  decoy 
D 


T  =  time  to  detect  an  HVS 
H 


with  odf's 


f  =  p  e 
H  hH 


V 


V 


t  *  0 


t  s  o 


Also  in  discussing  transition  probabilities,  reference  was  made  to  rate 

the  rate  at  which  an  HVS  encounters  decoys.  These  are  three  entities 

with  the  common  feature  of  respresenting  Poisson  processes,  that  is,  the 

time  between  detections  (or  encounters)  has  the  exponential  distribution. 

The  rates  corresponding  to  these  processes  are  3  ,  P  ,  and  \  .  It  is 

H  D  c 

characteristic  of  the  Poisson  process  with  rate  \  that,  if  n  is  the  num¬ 
ber  of  occurrences  in  time  At,  then 


p(n=0 ;  At)  =  1  -  X  •  At  +  o(At) 


p(n=l;At)  =  X  •  At  +  o(At) 


We  assume  the  rate  \  is  dependent  on  r  and  cp  ,  where  r  and  cp  are  the 

i  J  i  j 

range  and  heading  angle  discussed  in  Subsection  3.2.  In  this  case 


p(n=0;At)  = 


/,  /J  p(n=o|r  ,<p  ;At)p(r  )p(cp  ) 
i  i  3  3 

EE  [l  -  X(i,.j)At  +  (KAt^pCr^pO?  ) 
i  J 

1  -  X(i, j)p(ri>p(tp  )  At  +  o(At) 

i  j 


That  is, 


?? 


X(i,.j)p(ri)p(cpj)  , 


where  \(i,  j)  is  the  rate  obtained  at  range  r^  and  heading  angle  cp^ . 

The  conditional  rate  Mi,  j)  is  obtained  from  results  contained  in  Ref.  6, 


explained  as  follows. 


Consider  a  general  target  with  speed  v  and  detection  range  R.  The 


submarine’s  speed  is  u  on  a  heading  cp  and  it  is  located  a  distance  r 

J  * 


from  the  center  of  A  .  It  is  shown  in  Ref.  6  that  the  rate  at  which 

o 


targets  are  encountered  at  bearing  0  off  the  submarine  heading  is  given 


X(P;i,j)  =  f(3;u,v)  , 


where  6^  is  the  target  density  with  submarine  at  r  and  for  v  2  u 


f(P;u,v)  =  v  cos  1  ~  cos  P^  cos  3  +  n/u  -  v2  cos2  p 


and  for  v  >  u 


I  -1  /  v  \  f 2  2  sT 

f(P;u,v)  =  (  v  cos  —  cos  P 1  +  vu  -  v  cos  p 


-1  /  u\  ^  O  ^  -1  /u' 

when  -  cos  - s  P  »  -  cos  I — 

\  v/  \v; 


or  when  cos 


en 


5  cos 


-1  /  u 


K) 


when  P  5  -  cos  1  or  p  ^  cos  1  (-  — ^ 


These  expressions  assume  an  unbounded  area.  The  effect  of  a  bounded  area 

such  as  we  are  dealing  with  is  to  set  X(P;j,  j)  =  0  when  the  point  on 

the  detection  circle  at  heading  P  is  outside  Aq,  since  target  density  6 

is  assumed  to  be  zero  outside  A  .  The  total  rate  \(i,  j)  for  a  submarine 

o 

on  heading  cp  at  range  r  is  obtained  by  integrating  \(p;i,  j)  with  re- 

J  A 

spect  to  P  over  the  applicable  range  of  P .  To  help  visualize  the  range 
of  P  an  example  case  is  shown  in  Figure  3.8.  In  this  example  \(i,  j) 


There  are  a  number  of  separate  geometrical  cases  that  require  individual 
attention  in  determining  \(i,  j).  The  logic  and  computations  to  handle 
these  cases  are  shown  in  flow  chart  form  in  Figure  3.9.  Note  that  in 
this  flow  chart,  to  conserve  space,  we  denote 


>'•’  (i,  j)p(rj)p(cpj) 


i  J 


where  X;(i,  j)  is  obtained  from  X(i,  j)  with  n 


In  terms  of  earlier 


notation  we  have  then 


■  V« 

®D  =  Vd 

Xc  ’  "hY> 


where  X  ,  X^,  and  X^  are,  respectively,  the  rate  of  submarine  encounter 

of  a  single  HVS,  the  rate  of  submarine  encounter  of  a  single  decoy,  and 

the  rate  of  single  HVS  encounter  with  a  single  decoy.  Components  of 

p (cp  )  are  also  included  in  Figure  3.9.  We  take  p(cp.)  as 
J  J 


where  1 .  is  the  distance  from  the  submarine  along  its  heading  cp .  to  the 
J  J 

closer  of  two  boundaries:  (1)  the  detection  circle  R  centered  at  the 

submarine,  and  (2)  the  boundary  of  the  operating  area  A^.  Thus,  p(cp^) 

is  an  aa  hoc  measure  of  the  likelihood  of  the  submarine's  being  on  a 

particular  heading  cp,.  Basically,  this  measure  discounts  the  probability 

3 

of  the  submarine's  being  on  a  heading  directed  out  of  A  when  he  is 

o 

close  to  the  boundary.  The  probability  p(r_^)  is  determined  as  before, 


i.e.,  p(r^)  =  w(r^) . 


Conditioning  on  r.  and  cp .  and  computing  X  as  a  weighted  average 
were  introduced  to  account  for  the  facts  that:  (1)  we  are  dealing  with 
a  bounded  area,  and  (2)  when  the  submarine  is  in  the  search  process  and 
has  no  contacts  the  target  density  in  the  unswept  area  is  affected  by 
the  size  of  the  swept  area.  It  is  interesting  to  compare  the  weighted 
average  X  with  the  unweighted  formulation,  call  it  X  ,  where 


X 


* 


n  R  n 

“2  f  f(p;u,v)dP 
o  *  o 


X*  is  the  formulation  used  in  the  earlier  work  (Ref.  3).  A  comparison 
* 

between  X  and  X  is  shown  in  Figure  3.10.  The  deviation  between  the 
two  is  most  extreme  for  larger  values  of  detection  range  R,  as  would  be 
expected.  Deviation  also  increases  with  target  speed.  An  interesting 
comparison  between  X.  and  X*  is  in  the  computation  of 


Consider  as  an  example  n^ 
200.  Table  3.1  shows  the 


n  X 
D  H 


12 


n  X  +  n  X 
H  H  D  D 


=  1,  R  =  50,  n  =  1 
,  H  ’  D 

comparison. 


and  10, 


and  R 

D 


100  and 


Table  3.1 

COMPARISON  OF  DETECTION  RATE  FORMULATIONS 


(P12(X*> 
Table  Entry  < 


R  =  100 

D 

R  =  200 

D 

%  =  1 

.6689 

.8011 

.6691 

.8650 

n  =  10 

D 

.9528 

.975 

.9528 

.984 

The  differences  indicated  by  this  table  range  from  small  to  nonexistent. 
It  seems  appropriate  that  in  future  work  a  simple  approximation  to  X 
based  on  Figure  3.8  should  be  adequate  for  most  situations. 
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FIGURE  3.10  DETECTION  RATE  VERSUS  DETECTION  RANGE 


4.  MODEL  SOLUTION  AND  COMPUTATION 
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Given  the  measure  of  effectiveness  (MOE)  and  the  details  of  model 
structure,  it  remains  to  show  how  the  MOE  is  related  to  the  model  and 
how  it  is  computed.  Recall  that  the  MCE  is  expressed  in  terms  of  the 
function  (p^(t),  the  probability  that  at  time  t  the  process  is  in  State  j 
given  that  at  time  zero  it  was  in  State  i.  Assuming  that  at  time  zero 
the  process  starts  in  the  entry  state,  we  easily  narrow  the  scope  to 

Further,  in  the  six-state  model  we  are  interested  in  miselassi- 
fication  of  a  decoy  because  the  submarine  is  removed  as  a  threat  so  here 

wo  focus  on  cp  (t)  and  ©  (t). 

19  110 

Sinee  the  six-  and  ten-state  models  entail  different  assumptions, 
they  have  different  methods  of  solution.  First,  the  six-state  model  is 
treated  and  then  the  ten-state. 


from  node  1  to  the  output  node  of  the  modified  flow  graph  is  s©  (s). 

i 

e  -1.^ 

Hence,  tp  (s)  -  s  t  (s).  This  transform  turns  out  to  be  a  ratio  of 
X  *j  Xo 


polynomials  which  can  easily  be  inverted  by  finding  the  roots  of  the 


denominator.  The  modified  process  flow  graph  is  shown  in  Figure  4.1. 


The  transmission  from  node  i  to  j  of  a  flow  graph  is  defined  as 


I  y*  -*k  Ak 

A  k^P  ij  ij  ' 
kGPij 


where  A  is  one  plus  the  sum  of  all  loop  products  of  all  loops  in  the 


graph,  P  .  ic  the  set  of  paths  leading  from  i  to  j,  t  .  is  the  path 
^  th 

transmission  of  the  k  path  leading  from  i  to  j,  and  A  .  is  equal  to 

ij 


one  plus  the  sum  of  the  loop  products  of  all  loops  that  share  no  node 


with  the  k  path  from  i  to  j.  In  our  case  t  .  is  a  function  of  s  so 

"J 


we  write  t  (s).  The  simple  loop  products  for  the  graph  in  Figure  4.1 

ij 


are:  -P,  -GF,  -GHK,  -DFQ,  -DHKQ.  They  are  negative  by  definition. 


The  multiple  loop  products  are  obtained  by  considering  all  possible 


combinations  of  simple  loops  that  have  no  nodes  in  common  and  multiply¬ 


ing  their  products.  The  multiple  loop  products  for  Figure  4.1  are: 


FGP,  GHKP.  Hence,  we  have 


A  =  1  ~  [P  +  GF  +  GHK  +  DFG  +  DHKQ]  +  [FGP  +  GHKP] 


There  are  seven  paths  from  node  1  to  node  3.  They  can  be  identified  by 


their  transmissions  t  =  B,  AC,  ADE,  AGM,  AGHN,  ADQM,  ADQHN.  The  sums 
k 

A  are  one  lor  all  paths  except  AC,  AGM,  and  AGHN  for  which  they  are 

X  o 


l-P  and  B  for  which  it  is  A.  Hence, 


~t  (s)  =  “  [BA  +  AC (l-P)  +  ADE  +  AGM(l-P)  +  AGHN(l-P)  +  ADQM  +  ADHQN] 

13  A 


=  B  +  -  [DE  +(C  +  GM  +  GHN) ( l-P)  +  DQfM+HN)] 

A 


1 


A  =  P12h*(s) 
B  "  Pl3h*(,( 

c  -  p23h*(») 
o  -  P26h*(t| 

E  “  P63h6*J* 
F  -  P42h*4(5) 

G  °  P24h2<Sl 


H  ‘  P4Bhi,S) 

K  «  P52h*(s) 
M  -  p43h*(s) 
N  -  P53h*($) 

P  "  p66h6<sl 
Q  “  P64h6*S* 


FIGURE  4.1  MODIFIED  PROCESS  FLOW  GRAPH  FOR  SIX  STATE  MODEL 


Each  of  the  letters  in  A  and  t  (s)  is  identified  with  some  p  h  (s) 

Jl  o  lj  1 

To  obtain  t  as  a  function  of  s  we  must  specify  the  forms  of  h^(s). 
i  =  1,  4,  5,  h^t)  is  of  the  form 


h^t)  =  o^e 


-V 


t  S  0 


which  implies 


h®(s) 


s  +  c* 


For 
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t  i  0 


For  i  =  2  and  6,  h_^(t)  is  of  the  form 


oyV  V1  -”iV 

Vt5  ■  ispirr1  e 


which  implies 


he(s)  =  - - - 

i  s  +  no  | 

1  1  u 


To  simplify  algebra  at.u  subscripting  the  following  identifications  are 


O'  =a=n\+n\ 

1  H  H  D  D 

“2  “  b  *  ““d 

"«  - d  -  m 

"s’1-  Yh  +  Vi> 

as  -  c  ’  %  • 

£  “1  -As 

Now  cp  (s)  =  s  t’  (s)  can  be  written  in  terms  of  a  ratio  F(s)/Q(s) 

i.O  id 


of  polynomials  P(s)  and  Q(s),  specifically, 


‘O--0 


B  A  P(s) 

S  S  Q(S> 


where 


Np(n) 

7'  i-I 

,  V 


Nn(n) 


=  t  v 


In  tne  case  n  =  1,  we  find  that  N  (n)  =  4.  N  (n)  =  5,  and 

P  '  Q 
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V~.^55E? 


^KW^y-.’* v  ?  ~**'Z  "~v.  >  vK'^--'  n  j.  "  ..*y ~ , 


and 


P4  =  bp23 


=  '“'Wea  +  b(d+£)p23  +  bdP24p43 


b°<dt,>Cp26P63  +  a-p66,p23] 


+  bcd[p  p  (l-p  )  +  p  p  p  1 
24  43  66  *26*64P43J 


+  bcjf  (p„„+p„  p  +p  p  p  ) 
23  *24*43  *24*45*53 


bCdlCp26P63  +  P26P64(P43+P45P53) 


+  (1“P__)(P„  +P  P  +P  P  P  )] 
66  *23  24*43  P24P45P53/J 


%  =  1 


=  b  +  c(l-p  )  +  d  +  f 
66 


<lo  = 


=  (b+drf)c(l-pgg)  +  bd(l-p24p42)  +  (b+d)f 


bcd(l-p  -p  p  +P  P  W  -p  p  p  ) 
66  24  42  24  42*66  P26P42P64' 


+  (b+d)cf (1-cp  ) 
66 


+  bdf(l-p  p  -p  p  p  ) 

24  42  P24P45P52^ 


bcdf  (1-p  -p  p  -p  p  p 

66  *24*42  *24*45P52 


P26P42P64~P45P52P26P64+P24P42P66 


+  P24P45P52P66* 


For  cases  with  n  >  1  the  algebra  involved  in  obtaining  P(s)  and  Q(s)  is 
extremely  cumbersome.  For  this  reason  a  computer  program  utilizing  the 
IBM  FORMAC  algebraic  language  has  been  developed  to  handle  the  algebraic 
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manipulations.  FORMAC  is  an  extension  of  PL/1  and  provides  for  the 
symbolic  manipulation  of  mathematical  expressions,  e.g.,  the  expression 
SIN(X)  can  be  differentiated,  resulting  in  the  expression  COS(X),  Ex¬ 
pressions  can  be  differentiated,  evaluated,  replaced,  compared,  and 
parsed.  The  language  is  explained  in  detail  in  Ref.  7  and  the  program 
developed  for  this  project  is  listed  in  Appendix  D.  The  results  of  the 
program  in  terms  of  P(s)  and  Q(s)  for  n  =  2,  3,  and  10  are  listed  in 
Appendix  E. 

e 

Returning  to  the  expression  for  cp  (s), 

X  o 

e  BA  P(s) 

CD  (s)  =  ~  + - 

V13v  s  s  Q(s) 


and  substituting  for  A  and  B,  we  obtain 


me  ,  .  1  P(s) 

13  S  **13  s(s+a)  +  P123  s(s+a)Q(s) 


The  problem  now  is  to  invert  the  exponential  transforms.  For  the  first 
term  on  the  right  we  have  simply  from  a  table  of  transforms 


13  s(s+a) 


-at 

P13a-e  ) 


The  second  term  on  the  right  is  more  involved.  Let 


T(s)  =  s(s+a)Q(s) 


W(s)  =  P(s)/T(s) 


We  can  numerically  find  the  roots  of  Q(s)  since  we  know  the  coefficients; 

call  these  roots  r  ,  ...,  r^,  some  of  which  may  be  complex.  These  are 

also  roots  of  T(s).  T(s)  in  addition  has  roots  r  ,  =  -  a  and  r  =0. 

N+l  N+2 

Now  it  is  known  (e.g.,  Ref.  8)  from  algebra  that  W(s)  can  be  decomposed 
into  partial  fractions 


-- 


P 

% 


1 


N+2 


W(s)  = 


s  -  r. 


where 


P(r  )/T'  <r  ) 
k  k 


P(r  ) 
k 


r  (r,  +a)Q'  (r  )  +  r  Q(r  )  +  (r  +a)Q(r,  ) 


k  k 


k  k 


Then  inverting  \V(s)  we  get 


^  rkt 
w(t)  =  y  Z  e 
k 


So 


P  & 

“l  O 


P(s) 


12  s(s+a)Q(s) 


N+2  r  t 

pi2'  J2  v k 


Putting  all  the  terms  together  we  have  for  the  six-state  model 


N+2  r  t 

V’  ’  P13(1  '  6  >  +  Pl2a  4,  V  “ 


4.2  The  Ten-State  Model 


The  ten-state  model  is  restricted  to  exponential  type  pdf's,  so 
more  efficient  methods  of  solution  are  available  for  it.  The  exponential 
pdf’s  insure  that  the  ten-state  model  is  a  Markov  process.  We  are  in¬ 


terested  in  cp  (t)  and  cp  (t),  or  equivalently  cp_(t)  and  cr  (t)  in  the 
xy  i  iu  y  xu 


vector  cp(t)  defined  by 


<p(t)  =  cp(0)f(t)  , 
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where  $(t)  =  [to  .(t)]  and  9(0)  -  (1,  0,  0,  ...,  0).  Howard  (Ref.  4)  shows 
At  ^ 

that  $(t)  =  e  (matrix  exponential)  where  A  =  A(P  -  I)  with  A  a  diagonal 
matrix  with  elements  corresponding  to  the  rates  of  the  exponential  hold¬ 


ing  times.  Specifically,. 


where 


A  =  diag(X^) 


p  +  p 
H  D 


X  =  n  X  +  (n  -1)X  +  X 
6  H  H  D  D  R 


X7  -  WVd*1. 


X  =  n  X  +  n  X 
8  H  H  D  D 


X9  =  ° 


X  „  =  0 

10 


An  iterative  scheme  is  used  for  calculating  the  vector  9(t)  Given 
At 

9(t)  =  9(0)e  we  can  write 

...  A(t+At)  ...  At  AAt  .  .  A* At 

cp(t+At)  =  cp(0)e  =  cp(U)e  e  -  cp(t)e 


Thus,  to  compute  cp(t)  at  points  t  that  are  At  apart  we  need  only  compute 

AAt  Z 

one  matrix  exponential,  e  .  To  compute  the  matrix  exponential  e  we 

use  the  Pade  approximation 

Z  -1 

6  Z1  Z2  * 


where 
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5.  LOGICAL  FLOW 


This  section  summarizes  the  preceding  sections  on  model  structure 
and  geometrical  analysis  in  terms  of  logical  flow  charts.  Figure  5.1 
gives  a  summary  flow  chart  for  the  six-state  model  and  Figure  5.2  for 
the  ten-state  model.  These  charts  show  the  relationships  among  the  in¬ 
puts  required,  preliminary  calculations,  the  major  loop  on  number  of 
decoys,  transition  probability  and  rate  determination,  partial  fraction 
expansion  or  matrix  exponentiation,  and  calculation  of  trapping  proba¬ 
bilities. 
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Figure  6.1  shows  two  cases  of  model  comparison  involving  decoys; 
the  variable  parameter  is  decoy  detection  range.  Submarine  speed  is  10 
knots;  the  other  parameter  values  are  shown  below.  Case  1  is  a  single 
CV  in  a  nonconvergence  zone  with  a  larger  number  of  lower  fidelity 
devices;  Case  2  is  an  integrated  CTG  in  a  convergence  zone  with  a  sr^ller 
number  of  higher  fidelity  devices.  In  this  comparison  there  is  consider¬ 
able  difference  between  the  models.  The  maximum  difference  in  Case  1 
is  0.20  (difference  between  0.19  and  0.39);  the  maximum  difference  in 
Case  2  is  0.25  (difference  between  0.32  and  0.57).  It  will  be  shown 
below  that  the  difference  is  due  to  the  zero  speed  of  the  decoys,  which 
implies  a  low  rate  of  dissipation  of  swept  area. 


Case  1  Case  2 


CV  in  NCZ  CTG  in  CZ 


The  vertical  broken  lines  in  Figure  6.1  indicate  decoy  ranges  which 
are  subject  to  sensitivity  analysis  on  number  of  decoys  (Case  2  only) 
in  Figure  6.2.  The  value  =  70  nmi  is  the  break  point  between  the 
two  models  in  Case  2  of  Figure  6.1,  and  the  differ erne  net  uun  models 
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is  insensitive  to  decoy  numbers  in  Figure  6.2,  On  the  other  hand, 

R^  =  160  nmi  is  a  point  of  wide  disparity  in  Figure  6.1,  and  Figure  6.2 
indicates  this  disparity  is  a  function  of  number  of  decoys,  with  disparity 
maximum  at  four  decoys  and  decreasing  with  fewer  or  more  decoys. 

Thus  far  it  is  apparent  that  geometry  representation  can  have  a 
large  effect,  but  the  mechanism  of  the  effect  is  not  clear.  Case  3 
shown  in  Figure  6.3  clarifies  this  mystery;  Case  3  is  a  high  speed  CV 
in  a  convergence  zone.  The  important  fact  illustrated  in  this  figure 
is  th  .t  the  large  difference  evident  between  the  new  and  old  models  at 
a  decoy  speed  of  zero  collapses  when  the  decoy  speed  is  increased  to 
10  knots. 


The  conclusions  are:  (1)  geometry  representation  can  have  an  ap¬ 
preciable  effect,  (2)  the  effect  depends  on  the  parameters  of  the  problem 
(e.g.,  decoy  detection  range,  speed,  capture  time,  and  number),  and  (3) 
the  effect  is  due  primarily  to  the  representation  of  the  dissipation  of 
swept  area. 


6.2  Probability  of  Misclassif ication  and  Information  Retention  Time 

Inclusion  of  misclassification  probability  and  retention  by  the 
submarine  of  information  on  the  last  decoy  classified  are  the  other  two 
major  modifications  in  addition  to  geometry.  Figure  6.4  illustrates 
the  highly  significant  effect  of  the  misclassification  probabilities 
p  and  p  .  However,  it  has  been  found  that  retention  time  of  inforrna- 
tion  has  virtually  no  effect  over  a  wide  range  of  situations  for  almost 
any  reasonable  retention  time  (e.g.,  0  to  1^0  hours). 
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7.  CONCLUSIONS  AND  RECOMMENDATIONS 


The  models  discussed  in  this  document  are  implemented  as  operational 
computer  programs  and  have  been  used  in  concept  evaluation  studies  re¬ 
ported  elsewhere  (Ref.  1).  Experience  in  the  development  and  use  of 
these  models  allows  some  general  conclusions  to  be  drawn  about  the 
methodology. 

Two  conclusions  are  easy  to  draw:  (1)  potential  for  target  mis- 
classification  is  a  highly  significant  part  of  the  problem,  and  (2) 
retention  of  information  on  the  last  target  classified  is  insignificant. 

Another  conclusion  is  that  FORMAC  is  an  interesting  and  powerful 
tool  for  analyzing  continuous  time  semi-Markov  process,  but  it  can  be¬ 
come  cumbersome  and  time  consuming  in  itself.  Its  use  is  best  restricted 
to  a  small  number  of  states  (certainly  no  more  than  six),  or  to  cases 
of  very  sparse  transition  matrices,  where  it  can  be  very  useful  in  de¬ 
riving  compact  closed  form  solutions.  However,  for  a  larger  number  of 
states  it  appears  a  discrete  time  approximation  is  the  better  approach 
where  consideration  of  non-Markovian  processes  is  desirable.  In  the 
work  described  here  it  seems  a  2-hour  discrete  time  step  would  be  the 
minimum  necessary  and  perhaps  longer,  say  5  hours,  would  be  acceptable. 

Very  significant  is  what  has  been  discovered  about  the  importance 
of  geometry.  Geometry  representation  can  have  an  appreciable  effect, 
but  the  effect  is  due  primarily  to  the  representation  of  the  dissipation 
of  swept  area.  The  following  simplified  approach  is  suggested  for  future 


work: 


(1)  Ignore  boundary  conditions 


(2)  Auopc  simple  modification  of  unbounded  area  detection 
rates  via  approximations  based  on  Figure  3.10  or  an 
extension  of  this  figure. 


Appendix  A 


DISTRIBUTION  OF  TIME  IN  THE  ENTRY  STATE 


The  entry  state  is  set  up  to  allow  consideration  of  nonstationary 
(time  varying)  encounter  rates  as  the  submarine  enters  the  operating 
area.  Let 


Np(t) 


=  encounter  rate  for  HVSs  at  time  t 


ND(t) 


=  encounter  rate  for  de<  oys  at  time  t 


V*> 


yt,  = 


N  (t)  +  N  (t) 

H  D 

ytN  (t) dT 
o  H 

y^N  _(T)dT 
o  D 

/t§(T)dT 

o 

time  to  detect  an  HVS  with  pdf  f^Ct) 

time  to  detect  a  decoy  with  pdf  f  (t) 

H 


Assume 


-\(t> 

fu(t)  =  N  (t) e  t  a  o 

n  H 


-y^ 

f^(t)  =  N^(t) e  t  a  0 


If  N  (t)  and  N  (t)  are  constants,  then  these  are  exponential  distribu- 
H  D 

tions.  The  holding  time  pdf  h^(t)  is  given  by 


7°\  /  So 
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TS 


h^(t)  =  §(t)e  t  S  0 

It  can  then  be  shown  that 

P12  =  P[TD  =  Bin(T0,TH)j 

CO 

-  f 

o 

To  obtain  N  (t)  and  N  (t)  wc  use  results  from  Ref.  6.  Consider 
H  D 

Figure  A.l,  where  the  submarine  is  heading  directly  in  toward  the  center 
of  the  operating  area.  Assuming  the  detection  circle  is  tangent  to  the 


1  -  e 


-V‘> 


N  (t)e 
n 


■v*> 


dt 


FIGURE  A.1  ENTRY  GEOMETRY  FOR  HVS  AT  TIME  t 

operating  area  at  time  zero,  the  distance  between  the  su' marine  ard  the 

center  of  the  operating  area  is  R  +  R  -  ut  at  time  t,  given  the  sub- 

II  o 

marine  speed  u.  The  angle  P  is  given  by 


’-r-.-^ji^-i^  ;v4*-v  ?■*<>&*&?  nw^sy 


P  (t)  =  arccos 

H 


2  2  2 
R„  +(R  +R  -ut)  -R 
H  Ho _ o 

2R„(R„+R  -ut) 

n  n  O 


Koopman  (Ref.  6)  shows  the  following.  Assuming  the  submarine  is  operating 
in  an  infinite  area,  the  encounter  rate  for  HVSs  experienced  in  the  bear¬ 
ing  region  (-P,  P) 


H 


26hrh  5 

— ~  f  f(P;u,VH)dp 
o 


where  6  is  target  density.  We  assume  this  same  formula  applies  approxi- 
H 

mately  to  the  situation  in  Figure  A.l,  so  that 

2Vh  /H<t> 

V*>  -  IT  J  f(6iU-Vde 

H  o 

where  A„  is  A  minus  the  overlap  of  the  detection  circle.  A  similar 
Ho 

expression  is  used  for  NQ(t).  There  are  a  number  of  logical  exceptions 
to  these  expressions,  but  they  are  not  discussed  here;  the  objective 
of  this  appendix  is  to  give  only  an  overview  of  the  approach  examined. 

The  next  step  is  approximate  h^(t).  The  approach  taken  is  to  fit 


siM  ~~  f  N 


(T)dT 


with  an  Erlang  distribution 


A 

/v  .  a  .  ,  -nAt 

H  (t;  X,n)  =  1  -  e 


n 


k=l 


a  k-1 

(Xnt) 
(k-1) I 


t  ^  0 


via  a  least  squares  technique  on  \  for  fixed  values  of  n,  and  then  select 
the  value  of  n  with  the  least  total  squared  deviation.  In  sixteen  cases 
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examined,  n  =  1  occurred  fifteen  times,  and  n  =  2  once.  The  n  =  2  case 
is  shown  in  Figure  A. 2.  Since  n  =  1,  the  exponential  distribution,  is 
so  heavily  predominant  and  the  procedure  so  time  consuming,  this  special 
consideration  of  entry  has  been  abandoned  and  entry  is  subsumed  under 
search. 


Appendix  B 


THE  ERLANG  DISTRIBUTION 


In  this  study  the  most  general  form  of  capture  time  distribution 
considered  is  the  Erlang  distribution.  Let  t  be  the  random  variable 
time,  n  >  0  a  positive  integer,  and  p  >  0  a  real  number.  The  Erlang 
density  and  cumulative  distribution  functions  are,  respectively, 

.  „  (lin)n  n-1  -nut 

g(t;n,n)  =  - — 7  t  e  ^  ,  t  S  0 

(n-1) : 


G(t;  n,  p.)  =  1  -  e  ■  ,  t  £  0 


with  the  convention  0!  *  1.  The  exponential  transform  is  given  by 


e ,  r  -st, 

g  (s;n,y)  =  ELe  J 


=  PH* 

[s+(in  J 


The  expectation,  variance,  and  second  moment  are 


Var(t)  = 


S (9 /  SI  Of^o i.  SS 

C  (2>\<wvK. 


By  varying  n  and  p  a  large  number  of  different  shapes  can  be  obtained, 
as  illustrated  in  Figure  B.l.  In  particular,  when  n  =  1  the  Erlang 
distribution  is  the  exponential  distribution,  and  when  n  -»  05  then 
Var(t)  -»  0,  which  implies  constant  capture  time.  Replacing  (in  by  p  in 
the  density  and  distribution  functions  yields  the  gamma  distribution. 


STANDARD  DEVIATION  (O  ) 


Appendix  C 


DISTRIBUTION  OF  RANGE 
FROM  CENTER  OF  OPERATING  AREA 

In  obtaining  expressions  for  P(r  |c  )  and  P(r  If  ),  a  pdf  k(r)  for 

XX  1  a 

the  distance  r  from  the  center  of  the  operating  area  A  was  used.  This 

o 

pdf  has  the  form 

2tt 

k(r)  =  —  r  0  ^  r  £  R  , 

A  o 

o 


where  R  is  the  radius  of  A  .  The  form  follows  from  the  assumption 
o  o 

that  position  is  uniformly  distributed  throughout  A^.  Let  (x,  y)  be 
the  rectangular  coordinates  of  position  measured  from  the  center  of  Aq. 
The  assumption  may  be  stated 


f(x,y) 


2  2  2 
if  x  +  y  SR 

o 

otherwise 


We  want  to  derive  k(r), 


the  pdf  cf 


/ 2  2 
r  =  V  x  +  y  . 


Let 


or 


u 
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x2  +  y2 
x2 


(C-l) 


(C-2) 
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where  the  sign  on  the  radical  depends  on  the  quadrant  of  A  (see  Figure 

o 


C.l).  The  Jacobian  of  the  transformation  (C-2)  in  quadrant 


A  A 

2  1 


A,  A, 

3  4 


FIGURE  C.l  THE  FOUR  QUADRANTS  OF  A_ 


A^  is  the  determinant 
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and  J 

ruw  -  wz 


Similarly,  J  =  -J  ,  J  =  J  ,  and  J  =  -J  ,  so  |  J  |  =  | j|  is  the  same 

6  X  o  I  t  X  jL 


over  all  quadrants.  The  transofrmation  (C-l)  takes  quadrant  A^  into  a 


region  in  the  uw-  plane.  This  mapping  is  illustrated  in  Figure  C.2 
and  is  described  by  the  limits: 


x:  0  *♦  R 


u:  0  -»  R 


y:  0  -♦  R 


The  quadrants  A  ,  A  ,  A  ,  also  map  onto  regions  B  ,  B  ,  and  B  ,  which  are 
2  3  4  2  3  4 


coincident  with  B^.  Consider  B^  the  image  of  A^;  we  can  use  the  theorem 


on  transformation  of  random  variables  (e.g..  Ref.  10)  to  write 


n&e*-.re&i 


FIGURE  C.2  THE  MAPPING  OF  A,  ONTO  B, 


It  is  clear,  however,  that  B  is  also  the  image  of  A  ,  A  ,  A  so  that 

X  Z  u  4 

considering  all  the  transformed  probability  density  we  have 


h(u,w)  =  —  —j==g  for  (u, w)  eB^ 

o  /uw-w 

Now  we  integrate  out  w  in  order  to  obtain  the  marginal  pdf  g(u)  of  u: 


u 

g(u)  =  /  h(u,w)dw 


=  TT 

T 


for  0  ^  u  ^  R 
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Appendix  D 


FORMAC  PROGRAM  LISTING 


This  appendix  presents  a  listing  of  the  program  written  in  the  IBM 
FORMAC  language  to  handle  the  algebraic  simplification  of  the  flow  graph 
transmission  for  the  six-state  model.  A  sample  output  is  also  given. 

The  listing  is  presented  in  order  to  give  the  interested  reader  who  is 
unfamiliar  with  FORMAC  an  idea  of  what  the  language  is  like  and  what  it 
can  do.  This  appendix  can  only  be  considered  as  an  accompaniment  to 
Ref.  7  (the  FORMAC  handbook);  the  intention  is  not  to  provide  an  intro¬ 
duction  to  FORMAC  but  rather  an  example.  Neither  can  the  program  be 
considered  efficient;  indeed,  many  improvements  are  obvious  to  someone 
with  only  limited  knowledge  of  the  language.  The  following  correspondences 
should  be  noted  in  reading  the  program  and  relating  it  to  the  problem: 


PROBLEM  NOTATION 


PROGRAM  NOTATION 


a  through  f 
A  through  F 
s  +  a 

(coefficient  ii  P(s)) 


(coefficient  in  Q(s)) 


A 

P(s) 

Q(s) 


N  (n)  +  1 
y 


N2,  N3 

AA  through  PP 
A  through  F 
SPA 

CUS(i-l) 

CDS(i-l) 

DEL,  DELTA,  DEIH 
UP1 
DWN1 
LI 


(\C\  (VwcL\oo 
iMo/wK 
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It  is  assumed  in  this  listing  that  n  =  1  and  that  b  =  c  (the  simplest 
case  considered — the  one  given  in  the  main  text) .  It  is  interesting  to 
note  that  FORMAC  is  imbedded  in  PL/1,  so  that  it  is  possible  to  incorporate 
the  numberical  computation  into  the  same  program  as  the  algebraic  manipu¬ 
lation.  This  was  not  done  on  this  project  because  an  IBM  computation 
facility  was  not  conveniently  accessible,  and  PL/1-F0RMAC  is  unique  to  IBM. 


PROGRAM  LISTING 


1  oft MAC  TRY  t  UT* ) 5 


INPUT  TL  FjRv^c  PREPROCESSOR 

°Ril('.2s  PbLCRDURF  OPTIOnSIMAIN); 

F?lR  VAC...I.!('T  IONS ; 

OPT  SET  ( l  IM?LENGTH=72  I : 

PUT  PAGE ; 

PUT  LISM*  PROJECT  1316-5  .. 

CCL  POWER  PIN  FlXEOfJl); 

CCL  HPN  Bill  F!XED(3H  ; 

UCL  HPD  P. !  N  Fixroui); 

LCTIFNCI  FT  I*  (  K  1)  *U  2)  I  ***  111  /$(  3)  **J  ( 1 1) ; 
Nl*l{ 

N2=l: 

N3=  1 ; 

N4=l  I 

N&=  1 ; 


PEGIN: 

LET  I  AX  =rTl"M ",A,SPA)5 
HX  =  FT("N?''fd,SPh>  ; 

CX  =FT ( "N2" » B- SPB) 5 
«;X  -rT(«N4'',0tSP0)  ; 

FX  =  F  TC"N5",F,SPF)  )  ; 

pR.lNr_<UT(Ax;bx;cx;Dx;FX)  ; 


LET  t  AA 

=  pi?*ax; 

lib  = 

pf 3*AX5 

CC  = 

P23*flX5 

01)  = 

!>26*BX 

EC 

=  Pb3*cx; 

FF  = 

F4Z*l)X5 

GG  = 

P2**HX; 

HH  = 

P45*DX 

KK 

=  P52MX5 

•IN  - 

P4340X5 

NN  = 

P5J*FX{ 

PP  * 

P6f *CX 

00 

=  P64*CX  I 

» 

P«  IM_JUT  (AA}BHsCC}l)iJ}EE  5FF  {OUJWHJKK  \';l-P  JOO;  I  ; 

LFTl  DELIA  =  1  -  (PP  ♦  C.GM  F  ♦  G(.*HH*KK  *  i!!)»;-F*JU  ♦  DD*HH1PKK*CkI 
♦  (FF*Go*PP  +  GG*HH*KK*PP) 5 
CLi  =  C  TDC v (DLL  1  A )  J5 
PMNI.i'UT  (DEI  TA50F.L); 

Auburn  ruA)  ? 

ATrwiZi-(  A’iui  ; 

LF  T  l  AdO=OEL  )  : 

Ai'OM  IZE(  C  EL  1  5 
INDEX*  1 ; 

Go  Tti  ST  25 

ST3  s 


LIT (GFL=AdG) 5 
ATl.Ml  Zu  (  AoGI  5 
p3INr_0UT(U£L) 

LC 1  (  ur.lR  =  RlPLACL  (GFLiSRA  ,  S  A ,  SEP ,  $*"N2"#R, 

SPCtS+"N3"',C«S  »0,s»*'N4',*r.,SPf  ,S  +  "N';"wh  I  )5 
ATI  MZM-.  ELR  )5 

LEKT2  =  (f'Hor.C  +  (LC  ♦  G(.*NN  +  Gi'-HI-HOA  )*  (  I-PP ) 
«■  DC  POO*  (  '•N  +  hhPNNII/CU; 

T?i  ■=  CoDtrtJ12)  ); 

ATuNU-lOFL)  5 
Li  T  (  4GG=  I  ?P  )  5 
A1CM1Z!  M2PI  5 
!NLf x=25 
Co  ?r-  st 25 

ST  4 : 

IBM  I2k-A'-0)  : 

4fr.«UZB  (  AOGI  5 
P'-p.;T_'IJl  (T2»)  5 

MT(  Tt'  =  RTPLACF  (  T  2  B  ,SI  A  ,  S+ "Nl  "*  A  f  SI' H  ,  S  ♦"M  "  *li , 
SI'L f  S ♦,”N3" *L  t  S'*0 » S ♦"N4"  *(•,  SP  F  * S  +  "N (  "*f  I  15 
AT!  *  I  ZE  (  TZc  I  5 

|Bt(i/P1  =  F  X  PAM'  (NOM  ( 1  2  I);  L'KN  -  OENl.R  (  1  V  15 
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PROGRAM  LISTING  (continued) 


Ml  *  HIuHPMaIuPUM  ); 
imijfl  *  l  XPAND(UP1*S)  )  5 
NU  =  INTtGEP (Mi)  ; 
oo  mo  U’  nu; 
no«n*1; 

LET  ecus  (  "N"  )  xCnEFMUPl ,  ) )  { 

pmM_uuT(cus("N"i  > ; 

ATOMIZUCUSCN"!  }; 

INI); 

A  TLPiZt ( UPl )  ; 

LFT  a’fcNMFXPAN.MlwNl;  LI  -MGl-PUh  (PWN1 ,5 ))  J 
LEHDUNI  »  eXPANr(U*«M*S  )); 

PH!\T_UOT(H); 
ati  Mi/Ki:w,Ni ; 

NO  *  INTf.otrt  (LI): 

!)»>  MO  To  NO; 

n«*n+1; 

LFT(CI)S(  "N“»  =C«JbhF  «  U«M ,  ) : 

P«  JiNl_UUT«COS«"N*M  ) ; 

A) r  vi  zt  (CoS ( "f." ) ) ; 

Fnp; 

ATlMIZf  (OmM  ); 

PUT  LIST  ('  FIJI)  'IF  kLKI  I'  )  5 
CU  TO  STC; 

ST2 : 

LF  T ( AT  =  EXPAND  ( MJM  (  Abl.  ))$ 

L»T  *  rXPAf.D  ( DE-SuP  (  AOG  )  » S 
/ 1 1  >  =  SPA; 

Zl?)  *  SPb; 

Hi)  -  SPC; 

Z(4)  =  spu; 

z(bi  *  s»f  ): 

[t:  i  =  l  To  «; 

ATr*izr(A!>;ns>s 
LrT(AS=o;rtb*o;  Y=z(Mi«n  ; 

LFTIHPNY  =  H  IlihPUW  ( AT  »  Y )  « 

ipnv  -  li.mpum  at  »y  ); 

HPI)Y  =  KIGHPOwl  ul  ,Y)  ; 

IPi'Y  *  L  'wPOe  (  AT  *  Y  I  >; 

HPN  -  IMcUEMMPNYI  ; 

L  Pi\  =  INTFGtMLPNYU 

hpo  =  r-  nor- (I'Pdy  ) ; 

L-,D  *  !NTCi,r.  "(LPOY); 
l  ‘  LPN<  =  .»  In'  N  Go  ’.i  ill; 

kh; 

I1’  i.pocm  lMr’>i  G'i  To  Ml; 

tLSl ; 

It  lP\>sl)>0  Ihr N  P.'VilP  -  IPU5 

iLS-  I'lVM  =  lpn; 

i  •  k  =  1  '  •  1  •  MP  ■. ; 

LSHA-j  =  AS  dll.  pr  (  AT  .Y  *Y**(  "K"  -  “H+.CP")  ); 

LU  (  AS=  L/.rA.'4ll(  AS  )  J  i 
i  in; 

i-.  t  *  LPu  re  upo; 

Ll  I  (  =  I'S  ♦  O'rH  I  •’!  ,)  *v»K"  |*Y<f*l  "K"  -  ''n^vr^«)  |; 

L*  i  (i  s=t  xi'anou  si); 
r  no  ; 

AM  VIZ’  <  A I  ;■*  T)  ; 

l.»  H  AT  =  As:  1>I  =  •;>»; 

S  N  : 
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PROGRAM  LISTING  (concluded) 


END; 

LET (  AdG  *  AI/RT); 

IF  INDFXal  I  HEN  CD  lu  ST  3;  ELSE  U!  TH  STa  ; 

SI  E  s 

FND; 


SAMPLE  OUTPUT 


FPF'JFCT  m»i-5  ....  F-1KMAL  TF'Y  *-UT 


AX 

X 

A  /  SPA 

HX 

s 

M  /  S»U 

ex 

= 

0  /  SPH 

l/X 

X 

n  /  spo 

FX 

X 

F  /  SPF 

AA 

* 

PI?  A  / 

SPA 

ti  B 

r 

P13  A  / 

SPA 

cc 

z 

P?3  B  / 

SPH 

on 

X 

P26  B  / 

SPO 

f  r 

X 

P63  3  f 

SPH 

FF 

= 

P42  0  / 

SPO 

or. 

= 

P24  0  / 

SPB 

HH 

= 

P45  P  / 

sro 

UK 

* 

F  P52  / 

SPF 

PM 

= 

0  P43  / 

SPO 

f"N 

z 

F  FM  / 

SPF 

rp 

= 

•<  PM  / 

SPl! 

00 

r 

F.  P  44  / 

SPH 

OUT  A  =  r  P'*4  P?4  f  r>  Pet-  Pt<  /  t  SPF  SPl  SPl  )  -  F  MS  P26  |.  b 


Pf>4  PR?  /  «  SPF  SIT  $P:J  )  -  p  P4S  P?4  n  MS?  /  |  SPF  b 


P  bt<  /  SPP  ♦  F‘?4  14 


tfo  /  t  sm  spo 


p?f  i;  5  r*  4  /  ( 


pi  sp'  )  -  »•?*  u'i?  .  ./  /  (  sp;  spt  j  ♦  i 

If  10  ->4 

-  rt>(  S”!)  ’24  F**t  .)  SPP  rt  <  -  v’f  SPl'  P4?  P?«.  0 


m  «•  <  in 


P  Pf)C-  ♦  l 


F*.*>  I’?-*  0  SPl 


let  I  42  .;  S^t!  U  Po<  ♦  ( 


P4S  P<-4  0  SPD  I)  PSP  ♦  (  SPl  r  SPl 


SAMPLE  OUTPUT  (continued) 


2  2  3  fc  11 

SPF  SPD  SPB  )  SPB  )  SPO  SPB  )  SFC  SPI)  I  SPI)  5P0  ♦  SPF  SPO  SPe 

3  b  II 
»  /  <  SPF  SPD  SPB  ) 


CEL  »  (  F  P45  P 24  0  B  P66  P52  -  F  P45  P?6  D  B  Po4  P52  -  F  P45  **24  D 


2  2 

SPb  B  “S2  -  SPF  SPI)  SPd  a  P66  ♦  SPF  P24  P42  D  B  Ph4  -  SPT  P4?  P26  U  rt 


P64  -  SPF  °24  P4?  0  SPH  rt  ♦  SPF  SPC  SPB  J  /  (  SPF  SPD  SPB  > 


2  2 

T2R  «  (  -  SPF  r 24  D  B  Pbfc  P43  ♦  SPF  P26  D  B  P64  P43  ♦  SPF  P24  D 

2  2 

SPB  R  P43  -  F  P4S  P24  i)  d  P6o  P52  ♦  F  P4S  P26  D  6  P64  P53  ♦  F  P45 


P24  C  SPH  B  P53  -  P23  SPF  SPO  B  P6fc  ♦  P23  SPF  SPD  SPB  B  ♦  SPF  SPD  P63 

2  2  2 
P26  B  )  /  (  F  P4S  PP4  J  B  P6o  P52  -  F  P45  P26  U  B  Ffti,  p*j2  -  r  F4S 


P24  r  SPE  B  P52  -  SFF  SPC  SPP  B  Pfcfe  ♦  SPF  P?4  P42  D  B  i>Fl  -  SPF  P42 

2  2 
P26  L  0  P44  -  JPF  P 24  P42  0  SPB  b  ♦  SPF  SPD  SPb  > 

2  2  2 
Cl'S  (0 )  *  -  F  P24  D  B  Hfcb  p<,3  ♦  F  P26  U  B  P(4  P43  ♦  F  P?4  D  B  P4l 

2  2  2 
-  F  F45  P24  C  B  P66  P53  ♦  F  F4$  P?b  D  I)  Pf  <•  P53  ♦  F  P4t>  P24  D  0 

2  2  2 
PS3  -  P23  f  l)  B  P6fc  t-  Pbj  P2o  on  t  p?2  r  ()  6 

?  2 

CUSIH  *  -  P24  I)  d  P06  P43  ♦  P2b  I)  B  Pfc4  P43  ♦  F  P<-4  D  U  P4’'  «■  P24 

2  2  2 

0  0  P43  ♦  F  »4b  p 24  0  t)  P53  -  P?3  D  H  R(  o  -  P23  F  P  P64  ♦  P23  F  I) 

?  2  2  2 
H  ♦  P(>3  P2f>  t.  rt  ♦  P?l  i:  b  ♦  F  H6(  P2b  P  ♦  l>23  F  B 


CUS  ( 2 1  *  P24  l)  13  P43  -  °2i  B  P6b  ♦  P2J  l)  H  ♦  P23  F  E3  ♦  ?63  P26  l'  ♦ 


F23  l> 


CHS  I  3  >  =  H2?  u 


»“wnw*asa» 


SAMPLE  OUTPUT  (concluded) 


iv 


1 


m 


I 


IT 


I 

£ 


i 


B' 


1 


L 


i 


i 


I 


m 


& 

% 


f 

IT 


LI  =  A 


»  2 

COS  1C )  =  F  PAfi  H24  L  rt  Pho  P?>2  -  f  u*.‘>  P2*  l)  a  P<><-  Pf>2  -  r  n4?  P?a  ,> 


2  2  2  2 
«  PS?  ♦  F  P2*  P4Z  0  t)  P60  -  F  C  r»  P6(.  -  I  Pa2  P2t  1.  b  Pf.A  -  F 


2 

2 

P24  P4?  f  H  ♦  F 

0  B 

? 

2 

cost  19  =  -  f  PAS 

P^4  0  tt 

ps;  - 

F  0 

1*  P66  -l)F 

P66 

♦ 

P?4 

p42 

0  F 

2 

£ 

2 

PS  6  -  Pi)  P  66  - 

P42  P?6 

o  a 

Pl>4 

-  F  P24  Pa?  0 

t)  ♦ 

2 

r  n 

a  •* 

OP 

2 

2 

P24  P4?  IB  ♦  F 

B 

2 

COS ( 2  )  *  -OH  P66  -FB 

P(>6  - 

a 

P66  -  P?4  P42 

o  a 

♦ 

2  o 

B  ♦ 

c  f  B  ♦ 

F  C  ♦  i< 

CUSClJ  *  ?  h  ♦  C-  ♦  F  -  rt  Pob 
CPS<4)  =  l 

H,[  .if  ^pruxT 


Appendix  E 


POLYNOMIAL  COEFFICIENTS 

This  appendix  contains  the  coefficients  of  the  polynomials  P(s) 
and  Q(s)  for  n  =  2,  3,  and  10.  The  coefficients  are  presented  in  FORTRAN 
format  due  to  the  fact  that  they  have  been  printed  directly  from  verified 
punched  cards  in  order  to  eliminate  typographical  errors.  To  assist 
interpretation  of  the  coefficients  the  following  correspondences  between 
model  terminology  and  FORTRAN  are  noted: 


FORTRAN 

* 

A,B,C,D,F 
A2,B2,C2,D2,F2 
PIJ  (e.g.,  P26) 
P(i) 

Q(i) 


The  P(i)  and  Q(i)  are  given  in  Tables  E.l  through  E.3. 


111 


pm  * 


P(4)  * 

P(5)  * 
0(1)  « 


F£S«2vC  ;7V-"i  **  4  i^s **&&%?&?&** S&. ' K.'Sj^'r-'AW'**  *3j7^’t’;t  Ji 


^s*ssr> 


Table  E,1 


COEFFICIENTS  FOR  n  =  2 


l6*F*0#P26*C2*tf2*P53*P4S»P64 

♦  16-'F*D#P26*C2*82#P43<*P64 
-16«F»0#P24»C2*B2#P66»P53*P45 

♦  16»F*l)*P24«C2*B2*P53»P45 
-i6*F»D#P24*C2*92*P66»P43 
♦16*F»0*P24*C2*82*P43 
•16*F*0*P23*C2*B2*P6b 
♦16*F»D*P23*C2*B2 
♦16«F*D#P63#P26»C2*B2 

16#0#P26#C2*B2*P43*P64 

♦16#F»D*P24*C*B2*P53#P45 

-16»0*P24*C2*B2*P66*P43 

♦16*F*D*P24*C»B2»P43 

♦  l6*D*Pi4#C*!*B2*P43 
-1 6#0*P23«C2»B2»P66 
-lb*F»P23*C2»B2»P66 
♦16»F*D«P23*C*d2 
♦16«O»P23»C2#02 

♦16*F»P23*C2*B2 

♦16«0*P63»P26*C2*B2 

♦16»F«P63»P26*C2»B2 

♦•F«0*P24*B2*P53*P45 

♦16*D*P24*C*B2*P43 

♦4*F*0*P24*B2»P43 

-16*P23*C2*82*P6b 

♦16*D*P23*C*B2 

♦16*F*P23*C*B2 

♦4*F*0«P23*B2 

♦16*P23*C2*B2 

♦16*P63*P26*C2»B2 

4#D*P24»B2*P43 

♦  16*P23<*C*82 
♦4*0*P23#B2 
♦4*F*P23*B2 

4*P23*82 

■»16*F*D»P26»C2*B2*P52»P45*P6* 
•16#F*D*  >2*P26#C2#B2*Pb4 

♦  l6*F»0#P2‘»#C2*B2#Pb6*P52#P45 
-lb#F#u*P24*C2*B2#PB2»P45 

♦  j6*F»0»P24»P«<!»C2*B2»P6b 
-16*F*0*C2»B2»PB6 
-16*F»D»P24«»P42*C2*B2 
♦16*F*0*C2»B2 
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i 


Table  E.l  (continued) 


Q(j!  ■  *i6*D«P42*P26*C2*82*P64 

•16*F*0*P24»C*B2*P52*P45 

•16*F*D«C2*MP66 

♦16*D*P24*P42*e2*B2*P66 

-16*D»C2*B2*P66 

•  1<**F*C2*B2*P66 

♦1**F*0*C2*B  i 

-1&*F*D*P24*P42*C*B2 

♦16*F*0*C*B2 

•lfc*D*P24»P42*C2*B2 

♦16*D*C2*B2 

♦16*F*C2*B2 

0(3)  ■  <.4*F*0*P24*B2*P52*P45 
«16*D*C2*9*P66 
•lfc*F*C2*B*P66 
•16*C2*B2»P66 
-4*F*P*C2*P66 
*l6*F*i)*C*B 
♦16*D*C2*B 
♦16*F*C2*B 
-16*0*P24*P42*C*82 
♦16#D*C*B2 
M6#F»C«B2 
-4«F*0*P24*P42*B2 
♦4*F*D*82 
♦i6*C2*B2 
**«F#q#C2 

0(4)  »  •16*C2*B*P66 

-4*D*C2*P6b 

-4*F*C2*P66 

♦16*0*C»B 

♦16*F*C*8 

♦4*F»D*8 

♦16*C2*B 

♦4«F*0*C 

♦16*C*82 

-4*D*P24*P42*B2 

♦4*0*82 

♦4*F*B2 

♦4*0*C2 

♦4*F*C2 
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Table  E.l  (concluded) 


Q(5)  s  -4*C2*P66 

♦♦*0*8 

♦4#D*C 

♦F*0 

♦4*82 

♦4*C2 

0(6)  > 

•»4*C 

♦  D 

♦  F 

0(7)  «  i 
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- -— --  -  -  ^ 1  ■■tim  mi  i 


P(l)  «  729*F*D*P26«*C3*B3*P53»P45»P64 

♦729*F*D*P26*C3»B3*P43*P64 
-729*F*0*P24*C3»B3#Pbb*P53«P4S 
♦729*F»D*P24*C3*B3*P53»P45 
•729#F*D#P24*C3*83#P66*P43 
♦729*F*D#P24*C3*B3*P43 
-729*F*0*P23*C3*B3*Pb6 
♦729*F*D*P23«C3*B3 
♦729*F»D*P63*P2f>*C3*33 

p (2)  *  729*0*P26*C3*B3*P434P64 

♦729*F*D#P24#C2*83*P53*P45 

-729*QftP24<*C3*B3»P66*P43 

♦729»F40*P24*C2*B3*P43 

♦729*0*P24*C3»BJ»P43 

-729*0*P23»C3*B3»P66 

•729*F*P23*C3*B3*P66 

♦729*F«D*P23«C2»03 

♦729#0*P23*C3*B3 

♦729*F*P23*C3*B3 

♦729*0*P63*P2b*C3*B3 

♦729*F*P63*P26*C3*B3 

P (3)  s  243#F#0»P24#C*B3^P53*P45 
♦243*F*0*P24*C«B3*P43 
♦729*0#P24#C2*B3#P43 
-729*P23*C3*B3*P66 
♦243*F*0*P23*C*B3 
♦  729*l)*P23*C2*B3 
♦729*F*P23*C2»B3 
♦729*P23*C3#B3 
♦729*P63*P26*C3*H3 

p (4)  »  27*F»D*P24»B3*P53«P45 

♦243*L)*P24*C*B3*P43 

♦27*F»U*P24403*P43 

♦24J#0»P23*C*B3 

♦2*3*F*P23»C»W3 

♦27*F»0*P23*B3 

♦729#P23#C2*B3 

P(5)  *  27*0#PH4*B3*P43 

♦243*P23#C«B3 

♦27*0*P23*B3 

♦27«F4R23*B3 

p (6)  *  27*P23*B3 
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Table  E.l  (continued) 


0(1)  *  -729#F*l)*P26*C3*B3*P52*P45*P64 
-729*F*0*P42*P26*C3*B3*P&4 

♦  729*F#l)*P24*C3*B3*P66*P52«P45 
-729*F*U*P24*C3*B3*P52*P45 
♦729*F»D*P24*P42*C3*B3*P66 
-729*F*0*C3*B3*P66 
-729*F«'B*P24#P42*C3*B3 
♦729*F*0*C3*B3 

0(2)  *  -729*l)*P42*P26*C3«B3*P64 

-729*F»0*P24*C2»B3^P52*P45 

-729*F*D*C3*B2*P66 

♦729*0*P24«P42*C3*B3*P66 

-729*D*C3#B3*P66 

-729*F«C3*b3*P6ft 

♦729*F*D*C3*B2 

•729*F*0#P24*P42*C2*B3 

♦729*F#0#C2*H3 

-729#0*P24*P42*C3*B3 

♦729*D*C3*83 

♦729*F*C3*B3 

0(3)  »  -243*F#U*P24*C»B3*P52»P45 
-243*F*0#C3*B«P6b 
-729»0*C3*B«!*P66 
-729*F*C3*B2*P6fc 
•729*C3*B3*P66 
♦243*F*D*C3*B 
♦729*F»D*C2*B2 
♦729*D*C3*B2 
♦729#F»C34H2 
•243*F*D*P24*P42*C*B3 

♦  2-»3*F*D*C*B3 
-729*0«P24*P42*C2*B3 
♦729*D»C2*B3 
♦729*F*C2*B3 
♦729*C3*B3 


116 


a*-- 


*■*>*'■ 


'gH&Sp) 


Table  E.2  (continued) 


q<4)  ■  -27«F*0*P24»B3*P52*P45 
-243*0*C3*B*P66 
-243*F*C3»B*P66 
-729*C3«B2»P66 
-27*F#0*C3*P66 
♦24J*F»0*C2*B 
♦243*0*C3*B 
♦243*F»C3»B 
♦243*F*0*C*B2 
♦729*0*C2*B2 
♦729*F*C2*B2 
♦729*C3*S2 

•243*0*P24*P42*C*B3 

♦  243*l)*C*B3 

♦2*3*F*C*B3 

-27*F»0*P24*P42*B3 

♦27#F*0«33 

♦729*C2*B3 

♦27*F*D*C3 


U(5)  * 


-243*C3*B*P6b 

-27#D»C3*Pt»6 

•27«F*C3*P66 

♦81*F*0*C*B 

♦  243*0*C2»t» 

♦243*F*C2»B 

♦243*C3*B 

♦243*0*C*B2 

♦2*3*F*C*B2 

♦27*F»0*82 

♦729#C2#B2 

♦243»C*B3 

-27*0»P24«P42*B3 

♦27*0»B3 

♦27*F*B3 

♦27«F*0*C2 

♦27*0»C3 

♦27#F*C3 


Table  E.2  (coneludod) 


0(6) 


-27*C3*P66 

♦Hi#0*C#B 

♦81*F*C*B 

♦9*F*D»8 

♦2*3*C2*B 

♦9*F*D*C 

♦2*3*C*B2 

♦27*0*82 

♦27*F*B2 

♦  27*I)*C2 

♦27*F*C2 

♦27*83 

♦27*C3 


0(7) 


8l*C*B 

♦9*0*8 

♦9*F*8 

♦9*0*C 

♦9*F*C 

♦F*0 

♦27*82 

♦27*C2 


0(8) 


9*0 


♦9*C 
♦  0 
♦F 


Q(9) 
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Table  E.3 


COEFFICIENTS  FOR  n  =  10 


P(l)  *  ( -F*P24«U*B20*P66»P43 

♦F*P26*0*820*P64*P43 

♦F*P24*0*820*P43 

«F*P45*P24<»D*B20»P66*P53 

♦F*P45*P26*D*B20*P64»P53 

♦F*P45*P24*D*B20*P53 

-P23#F*0*820*P66 

♦F*P63#P26*0*820 

♦P^3*F*D*B20) •1E*20 

P(2)  s  (.P24*0«B20*P66*P43 

♦P26«0*820*P64*P43 

♦F*P24*D*B19*P43 

«pC4*0*B2o*P43 

♦F*P45*P24*D*B19»P53 

"P23*0#820#P66 

-P23*F*B20*P66 

♦P23*F*0#B19 

♦P63«P26*0«82o 

♦P23*D»B20 

♦F*P63*P26*B2f> 

♦  P23*F*B20>  *1E*20 

P(3)  *  (.45*F*P24«0*B18*P43 

♦P24*D«B19*P43 
♦ , 45*F*P45*P24*0*B i8*P53 
-P23*B20*P66 
♦,45*P23*F*D»B18 
♦P23*D*B19 
♦P23*F*B19 
♦P63»P26*820 
♦P23*B20)*1E*20 

P(4)  «  (,12*F*P24»0*B17*P43 

♦,45*P24*D*B18*P43 

♦,12*F*P45»P24*D*B17*PS3 

♦ » 1 2*P23*F#0#B1 7 

♦,45»P23*D*B1B 

♦.45*P23*F*B18 

♦P23»B1V)*1E*20 

P<5)  «  ( •21*F*P24*U#B16#P43 

♦1.2«P24«D*B17»P43 

♦.C1*F»P45»P24<»0*B16*P53 

♦  ,21<*P23*F<*0*B16 
♦l*2»P23*D*Bi7 
♦1.2*P23*F*B17 
♦4»5*P23*Bi8) *1E*19 
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Table  E.3  (continued) 


P(6>  ■  ( #252*F*P24*D*B15#P43 

♦2« 1*P24*D*B16*P43 

♦  ,252«F»P45«>P24*0*B15*PS3 

♦  ,252*P23<»F<»D*815 
♦2.1*P23*D*8H 
♦2.1*P23*F*B16 
♦12,*P23*B17)*1E*18 

P(7)  *  U21«F*P24»D*B14*P43 

♦2.52«P24*D»Bl5*P43 

♦.2l«F*P45*P24*0«Bl4oP53 

♦  «2l#P23#F*f)#814 
♦2.52*P23»D»B15 

♦  2.52»P23«*F*B15 

♦21 •*P23#B16) * 1 E ♦ 1 7 

P (8)  *  <.12*F*P24*l)»B13*P43 

♦2« 1#P24*D*B14*P43 
♦,12»F*P4S*P24*D*B13*P53 
♦• 12*P23#F*D*B13 
♦2« 1*P23#D#B1 4 
♦2« 1*P23#F#B14 
♦25.2*P23*Bl5) *l£*16 

P (9)  «  J,45*F*P24*D*B12*P43 

♦12,»P24*0*B13*P43 

♦,45*F*P45*P24*D*812*P53 

♦,45*P23«F*n«B12 

♦l2.*P23*0»Bi3 

♦12.*P23*F«*B13 

♦210.*P23*B14)*1E*14 

P ( 1 0 >  *  IF*P24*D*BU#P43 

♦45«*P24*D*B12#P43 
♦F*P45#P24#D*8l 1#P53 
♦P23#F*D*B11 
♦45«*P23*D*B12 
♦45»*P23#F#B12 
♦ 12OO»#P23*0 13) * 1 E ♦ 1 2 

P(ll)  *  ( F*P24*D*H 1 0#P43 

♦H)0.*P24»U*BU#P43 

*F*P45*P24*D*B10'*P53 

♦P23»F»D»B10 

♦  1 00«*P23#D#B1 1 
♦100.*P23*F*B11 

♦  4ijOO»#P23*B12>  *1E*  10 
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Table  E.3  (continued) 


P(12)  •  (P24*D»B10*P43 

♦P23*D#B10 

♦P23*F*8l0 

♦100.*P23*BH)»1E*10 

PU3)  ■  P23*B10*1E*10 

0(1)  *  (F«P45»P24*0*820*P66*P52 

-F*P45»P26*0*B20*PB4*P52 

-F*P45*P24*0*B20fP52 

♦F*P24*P42*0*B20*P66 

-F»0»B20*P66 

•F*P42#P26*D*B20#PB4 

-F*P24*P42*0*B20 

♦F*0*B20)*1E*20 

0(2)  «  (-F*P45*P24*0*B19»P52 

•F*0*B19*P66 
♦P24*P42«0*B20*P66 
-D*B20*P66 
-F*B20*P66 
-P42»P26»0*B2o*P6* 
•F*P24*P42«0*B19 
♦2«*F*0*819 
.p24»P42«0*B20 
♦0*B20 

♦F*820>*1E*20 

0(3)  ■  («,45*F»P45*P24#0*B1B*P52 

-,45*F»D*B18»P66 
•0*B19*P66 
-F*B19*P46 
-BC0#P6B 

-,45»F»P24*P42*0*B18 

♦1«9*F*0#B18 

»P24*P42*D*B19 

♦2.<*0*B19 

♦2«*F*8l9 

♦B20) *1E*20 


w 

b 


t 


q ( 4 )  «  (-.12*F*P45*P2**0*B17*P52 

«, 12*F*D*817*Pbb 
45*0*8 18*P66 
-,45*F*B18*P66 
-B19*P66 

-,12*F*P24*P42*0*B17 
♦ 1 « 14*F*D*Bl7 
-,45*P24*P42*D*Bl8 
♦1.9*D*B18 
♦1 »9*F*B18 
♦2»*B19) *1E*20 

q ( 5 )  s  («.21*F«P45*P2**0*Blb*P52 
-,21*F*D*B16*P66 
-1.2*0*B17*Pb6 
-1.2*F*B17*P66 
•4«5*B18#P66 
“,2i*F*P24*P*2*D*B16 
♦4.845*F*0*B16 
-1.2*P24*P*2*D*817 
♦ 1 1 «4*0*B1 7 
♦ll.**F*Bl7 
♦19,*B18)*1E*19 

0(6)  *  (-,252*F*P*5*P24*D*B15*P52 

-,252*F*0*B15*P66 
-2.1*0*B16*P66 
-2,1*F*B16*P66 
-l2,*B17*P66 
-,252*F*P24»P42*D*815 
♦15.5o^#P*0*B15 
-2.1*P2**P42*D*B16 
♦48,45*0*816 
♦*8,45*F*B16 
♦1 14**817) *1E* 18 

Q(7)  «  (-,21*F*P45*P2**0*B1**P52 

-,21*F*D*814*P66 
-2.52*0*B15*P66 
•2,52*F*Blb*P66 
-2l,*816*P6b 
-,21*F*P24*P42*0*B14 
♦38,76*F*0*814 
-2,52*P2*#P*2#0*B15 
♦l55,0**0*Bl5 

♦155.0*#F*B15 

♦484»5*Blb)*iE*i7 


t 


Table  E.3  (continued) 


0(8)  •  <-.12*F*P45*P24*0*B13*P52 

-,12*F*D*B13*P66 
-2.1*D#B14*P66 
-2.1*F*B14*P66 
-25,2*Bl5*P66 
-,12*F*P24#P42*D*B13 
♦77,52*F*D*B13 
-2.1*P24*P42*Q*B14 
♦387,6*0*814 
♦387»6*F*814 
♦155o.**B15)*1E*16 

Q (9)  «  (..45*F«P45#P2**D*B12*P52 

-,45*F*D*B12*P66 
-12,*D*B13*P66 
-12.*F*B13*P66 
-2lO.*Bl**P66 
-,45*F*P24«P42*D*B12 
♦1259.7*F*D*B12 
-12,*P24*P«»2*D*B13 
♦7752,*D*B13 
♦7752,*F*B13 
♦38760. *B14)*1E*14 

0  (10)  •  <-F*P45*P24*0*81i*P52 
•F*0*B1 1*P66 
-45.*D*B12*P66 
-45,*F*B12*P66 
-120d.«B13*P66 
•F*P24*P42*D*B1 1 
♦ 16796. *F*D*Bll 
-45.*P24*P42*Q*B12 
♦125970. •0«B12 
♦125970. *F*B12 
♦ 775200. *B 13) *1E*12 

Q  ( 1 1 )  *  UF*P45*P24*0*B10*P52 
-F*D*B10*P66 
-1  o.#0*Bll*P66 
"1U0.*F*B11*P66 
-4500.*B12*P66 
■F*P24*P42*0*B10 
♦164756, *F*D*B10 
-100. *P24*P4f *0*811 
♦1679600.*0*B11 
♦1679600. *F*B11 
♦12597000.*B12) *lE*io 


AC&ftrife * ... 


Table  E.3  (continued) 


0(12)  «  (-D*810*P66 

«F*810*P66 
-100.*Bil*P66 
♦16796. *F*D*B9 
*P24*P42*D*B1 0 
♦184756, *U*B10 
♦184756. *F*BlQ 
♦1 679600 .*81 1 ) *1E*1 0 

0(13)  •  <-bl0*P66 

♦1259.7*F*D*B8 
♦16796. *0*B9 
♦16796. *F*89 
♦184756. *Bl0)*lE*10 

0(14)  *  ( ,7752*F*D*B7 

♦12.597*0*88 

♦12.597*F*B8 

♦  167.96*B9)*1EM2 

Q ( 15)  *  <.3876*F*0*B6 

♦7.752*0*87 

♦7.752*F*87 

♦  1<J5.97#aa)*lE*ll 

0(16)  «  (15.504»F*D*B5 

♦387.6*0*86 
♦387.6*F*06 
♦7752.*B7) *l£*8 

0(17)  *  (4,845*F*0*B4 

♦155.04*0*85 

♦155.04*F*B5 

♦3876.*p6)*lE*7 

0(18)  «  ( 1 14,*F*0*B3 

♦4845. *0*84 
♦4845.*F*B4 
♦15504U.*B5)*1E^4 

0(19)  •  19000. *F*0*B2 

♦1140000, *D*B3 
♦1140000. *F*B3 
♦48450000. *84 
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Table  E.3  (concluded) 


Q (20)  • 


0(21)  - 


0(22)  « 


200.*F#t>»9 

♦19uQ0«*D*B2 

♦19000.*F*B2 

♦ll*0000»*B3 

200. •0»8 

♦200.***tt 

♦F*D 

♦l9000.»B2 

200. *8 
♦0 
♦F 


0(23)  *  1. 


